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Abstract: One of the most surprising consequences of quantum mechanics is the entan- 
glement of two or more distant particles. In an entangled EPR two-particle system, the 
value of the momentum (position) for neither single subsystem is determined. However, if 
one of the subsystems is measured to have a certain momentum (position), the other sub- 
system is determined to have a unique corresponding value, despite the distance between 
them. This peculiar behavior of an entangled quantum system has surprisingly been ob- 
served experimentally in two-photon temporal and spatial correlation measurements, such 
as "ghost" interference and "ghost" imaging. This article addresses the fundamental con- 
cerns behind these experimental observations and to explore the nonclassical nature of 
two-photon superposition by emphasizing the physics of 2 ^ 1 + 1 . 



1 Introduction 

In quantum theory, a particle is allowed to exist in a set of orthogonal states simultaneously. 
A vivid picture of this concept might be Schrodinger's cat, where his cat is in a state of 
both alive and dead simultaneously. In mathematics, the concepts of "alive" and "dead" are 
expressed through the idea of orthogonality. In quantum mechanics, the superpositions of these 
orthogonal states are used to describe the physical reality of a quantum object. In this respect 
the superposition principle is indeed a mystery when compared with our everyday experience. 

In this article, we discuss another surprising consequence of quantum mechanics, namely 
that of quantum entanglement. Quantum entanglement involves a multi-particle system in 
a coherent superposition of orthogonal states. Here again Schrodinger's cat is a nice way of 
cartooning the strangeness of quantum entanglement. Now imagine two Schrodinger's cats 
propagating to separate distant locations. The two cats are nonclassical by means of the 
following two criteria: (1) each of the cats is in a state of alive and dead simultaneously; (2) 
the two must be observed to be both alive or both dead whenever we observe them, despite 
their separation. There would probably be no concern if our observations were based on a large 
number of alive-alive or dead-dead twin cats, pair by pair, with say a 50% chance to observe a 
dead-dead or alive-alive pair. However, we are talking about a single pair of cats with this single 
pair being in the state of alive-alive and dead-dead simultaneously, and, in addition each of the 
cats in the pair must be alive and dead simultaneously. The superposition of multi-particle 
states with these entangled properties represents a troubling concept to classical theory. These 
concerns derive not only from the fact that the superposition of multi-particle states has no 
classical counterpart, but also because it represents a nonlocal behavior which may never be 
understood classically. 

The concept of quantum entanglement started in 1935 [1]. Einstein, Podolsky and Rosen, 
suggested a gedankenexperiment and introduced an entangled two-particle system based on 
the superposition of two-particle wavefunctions. The EPR system is composed of two distant 
interaction-free particles which are characterized by the following wavefunction: 

*(x 1 ,x 2 ) = J d Pl dp 2 S(p 1 +p 2 )e i ^-^/ n e i ^ n = S(x 1 -x 2 -x ) (1) 

where e %Pl ^ Xl ~ x °^ '' h and e W2X2 ^ n are the eigenfunctions with eigenvalues p\—p and p 2 = —p of 
the momentum operators p\ and p 2 associated with particles 1 and 2, respectively. X\ and x 2 
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are the coordinate variables to describe the positions of particles 1 and 2, respectively; and xq 
is a constant. The EPR state is very peculiar. Although there is no interaction between the 
two distant particles, the two-particle superposition cannot be factorized into a product of two 
individual superpositions of two particles. Remarkably, quantum theory allow for such states. 
What can we learn from the EPR state of Eq. ([1])? 

(1) In coordinate representation, the wavefunction is a delta function 5(xi — X2 — x§). 
The two particles are separated in space with a constant value of x\ — xi = Xq, although the 
coordinates x\ and x-i of the two particles are both unspecified. 

(2) The delta wavefunction S(x% — X2 — xq) is the result of the superposition of plane 
wavefunctions for free particle one, e Wl ^ Xl " x °^ h , and free particle two, e lp2X2 ^ h , with a particular 
distribution 5(p± + P2)- It is S(j>i + P2) that made the superposition special. Although the 
momentum of particle one and particle two may take on any values, the delta function restricts 
the superposition to only those terms in which the total momentum of the system takes a 
constant value of zero. 

Now, we transfer the wavefunction from coordinate representation to momentum represen- 
tation: 



*(pi, P2) = 7TT / dx ^ S ( x i ~ x * ~ *o) e -*i(«i-xo)/ft e -fcW* = 6(pi + pa). (2) 



What can we learn from the EPR state of Eq. @? 

(1) In momentum representation, the wavefunction is a delta function S(pi + P2)- The 
total momentum of the two-particle system takes a constant value of pi +P2 = 0, although the 
momenta p\ and P2 are both unspecified. 

(2) The delta wavefunction 5{jpx + P2) is the result of the superposition of plane wavefunc- 
tions for free particle one, e -' l Pi{xi-x )/h ^ an( j f ree p ar ticle two, e~ lp2X2 ^ h , with a particular 
distribution 5(xi —X2 — xq). It is S(x\ —X2 — xq) that made the superposition special. Although 
the coordinates of particle one and particle two may take on any values, the delta function 
restricts the superposition to only those terms in which x\ — X2 is a constant value of xo- 

In an EPR system, the value of the momentum (position) for neither single subsystem 
is determined. However, if one of the subsystems is measured to be at a certain momentum 
(position), the other one is determined with a unique corresponding value, despite the distance 
between them. An idealized EPR state of a two-particle system is therefore characterized by 
A(pi +P2) = and A(xi — x%) — simultaneously, even if the momentum and position of 
each individual free particle are completely undefined, i.e., Apj ~ 00 and Axj ~ 00, j = 1,2. 
In other words, each of the subsystems may have completely random values or all possible 
values of momentum and position in the course of their motion, but the correlations of the two 
subsystems are determined with certainty whenever a joint measurement is performed. 

The EPR states of Eq. JTJ) and Eq. ([2]) are simply the results of the quantum mechanical 
superposition of two-particle states. The physics behind EPR states is far beyond the acceptable 
limit of Einstein. 

Does a free particle have a defined momentum and position in the state of Eq. ([T]) and 
Eq. ([2]), regardless of whether we measure it or not? On one hand, the momentum and position 
of neither independent particle is specified and the superposition is taken over all possible values 
of the momentum and position. We may have to believe that the particles do not have any 
defined momentum and position, or have all possible values of momentum and position within 
the superposition, during the course of their motion. On the other hand, if the measured 
momentum (position) of one particle uniquely determines the momentum (position) of the other 
distant particle, it would be hard for anyone who believes no action-at-a-distance to imagine 
that the momenta (position) of the two particles are not predetermined with defined values 
before the measurement. EPR thus put us into a paradoxical situation. It seems reasonable for 
us to ask the same question that EPR had asked in 1935: "Can quantum-mechanical description 
of physical reality be considered complete?" [T] 
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In their 1935 article, Einstein, Podolsky and Rosen argued that the existence of the entan- 
gled two-particle state of Eq. |T]) and Eq. Q , a straightforward quantum mechanical superpo- 
sition of two-particle states, led to the violation of the uncertainty principle of quantum theory. 
To draw their conclusion, EPR started from the following criteria. 

Locality: there is no action-at-a-distance; 

Reality: "if, without in any way disturbing a system, we can predict with certainty the 
value of a physical quantity, then there exist an element of physical reality corresponding to 
this quantity." According to the delta wavefunctions, we can predict with certainty the result 
of measuring the momentum (position) of particle 1 by measuring the momentum (position) 
of particle 2, and the measurement of particle 2 cannot cause any disturbance to particle 1, if 
the measurements are space-like separated events. Thus, both the momentum and position of 
particle 1 must be elements of physical reality regardless of whether we measure it or not. This, 
however, is not allowed by quantum theory. Now consider: 

Completeness: "every element of the physical reality must have a counterpart in the com- 
plete theory." This led to the question as the title of their 1935 article: "Can Quantum- 
Mechanical Description of Physical Reality Be Considered Complete?" 

EPR's arguments were never appreciated by Copenhagen. Bohr criticized EPR's criterion of 
physical reality [2J: "it is too narrow" . However, it is perhaps not easy to find a wider criterion. 
A memorable quote from Wheeler, "No elementary quantum phenomenon is a phenomenon 
until it is a recorded phenomenon" , summarizes what Copenhagen has been trying to teach us 
[3]. By 1927, most physicists accepted the Copenhagen interpretation as the standard view of 
quantum formalism. Einstein, however, refused to compromise. As Pais recalled in his book, 
during a walk around 1950, Einstein suddenly stopped and "asked me if I really believed that 
the moon (pion) exists only if I look at it." [4j 



There has been arguments considering A(pi + £»2)A(xi — x 2 ) — a violation of the uncer- 
tainty principle. This argument is false. It is easy to find that pi + P2 and Xi — x 2 are not 
conjugate variables. As we know, non-conjugate variables correspond to commuting operators 
in quantum mechanics, if the corresponding operators exist0 To have A(f>i + p-i) = and 
A(xi — X2) = simultaneously, or to have A(pi + p2)A(xi — X2) = 0, is not a violation of the 
uncertainty principle. This point can easily be seen from the following two dimensional Fourier 
transforms: 



*(xi, x 2 ) = — ^ J dp 1 dp 2 S{pi+p 2 ) 



e ipi(xi—xo)/H gip2X2/h 



^ / d(pi+p 2 )8( P i+p 2 )e l ^ +P2){ < +X2)/2h I d( Pl -p 2 )/2e l ^- p ^<- x ^/ 2h 



= 1 x S(xi — X 2 — Xq) 
where x' = x\ — xq; 



ty(pi 7 p 2 ) = 7TT ( dxi dx 2 8{xi - x 2 - xq) e 
27m J 



—ipi (xi-Ko)/B e -ip2X 2 /fi 



= — f d(x' 1 +x 2 )e~ l{pi+P2){ < +X2)/2h [ d(x[ - x 2 )/2 6(x' 1 - x 2 ) e-' l{pi - p2) «- X2)/2h 
2-kTi J J 

= 5(pi +p 2 ) x 1. 

The Fourier conjugate variables are (xi +x 2 ) (pi +p 2 ) and (x\ —x 2 ) (pi — p 2 ). Although it 
is possible to have A(xi — x 2 ) ~ and A(pi +p 2 ) ~ simultaneously, the uncertainty relations 
must hold for the Fourier conjugates A(xi +x 2 )A(pi +p 2 ) > fr, and A(xi — x 2 )A(pi —p 2 ) > fr; 
with A(pi — p 2 ) ~ 00 and A(xi + x 2 ) ~ 00. 



It is possible that no quantum mechanical operator is associated with a measurable variable, such as time 
t. From this perspective, an uncertainty relation based on variables rather than operators is more general. 
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As a matter of fact, in their 1935 paper, Einstein-Podolsky-Rosen never questioned A(x± — 
X2) A(pi + P2) — as a violation of the uncertainty principle. The violation of the uncertainty 
principle was probably not Einstein's concern at all, although their 1935 paradox was based 
on the argument of the uncertainty principle. What really bothered Einstein so much? For 
all of his life, Einstein, a true believer of realism, never accepted that a particle does not have 
a defined momentum and position during its motion, but rather is specified by a probability 
amplitude of certain a momentum and position. "God does not play dice" was the most vivid 
criticism from Einstein to refuse the Schrodinger's cat. The entangled two-particle system 
was used as an example to clarify and to reinforce Einstein's realistic opinion. To Einstein, 
the acceptance of Schrodinger's cat perhaps means action-at-a-distance or an inconsistency 
between quantum mechanics and the theory of relativity, when dealing with the entangled EPR 
two-particle system. Let us follow Copenhagen to consider that each particle in an EPR pair 
has no defined momentum and position, or has all possible momentum and position within the 
superposition state, i.e., imagine Apj ^ 0, Axj ^ 0, j = 1,2, for each single-particle until the 
measurement. Assume the measurement devices are particle counting devices able to identify 
the position of each particle among an ensemble of particles. For each registration of a particle 
the measurement device records a value of its position. No one can predict what value is 
registered for each measurement; the best knowledge we may have is the probability to register 
that value. If we further assume no physical interaction between the two distant particles and 
believe no action-at-a-distance exist in nature, we would also believe that no matter how the 
two particles are created, the two registered values must be independent of each other. Thus, 
the value of X\ — X2 is unpredictable within the uncertainties of Ax\ and Ax 2 . The above 
statement is also valid for the momentum measurement. Therefore, after a set of measurements 
on a large number of particle pairs, the statistical uncertainty of the measurement on p\ + p 2 
and £1 — X2 must obey the following inequalities: 

A(pi +P2) = ^(Api) 2 + {A P2 ) 2 > Max(A Pl ,Ap 2 ) (3) 
A(xi - x 2 ) = ^{Axx) 2 + (A.t 2 ) 2 > Max(Axi,Ax 2 ). 

Eq. ^ is obviously true in statistics, especially when we are sure that no disturbance is possible 
between the two independent-local measurements. This condition can be easily realized by 
making the two measurement events space-like separated events. The classical inequality of 
Eq. (J3j> would not allow A(pi +P2) = and A(xi —x 2 ) — as required in the EPR state, unless 
Api = 0, Ap2 — 0, Ax\ — and Ax 2 = 0, simultaneously. Unfortunately, the assumption 
of Api — 0, Ap2 — 0, Axi = 0, Ax2 = cannot be true because it violates the uncertainty 
relations ApiAxi > h and AP2AX2 > h. 

In a non-perfect entangled system, the uncertainties of p\ +P2 and x\ — X2 may differ from 
zero. Nevertheless, the measurements may still satisfy the EPR inequalities [5]: 

A(p\ +P2) < min{Api, Ap 2 ) (4) 
A{x\ — X2) < min(Axi, Aa^)- 

The apparent contradiction between the classical inequality Eq. §3§ and the EPR inequality 
Eq. deeply troubled Einstein. While one sees the measurements of p\ +P2 and x\ —X2 of the 
two distant individual free particles satisfying Eq. (|4]), but believing Eq. ([3]), one might easily 
be trapped into concluding either there is a violation of the uncertainty principle or there exists 
action-at-a-distance. 

Is it possible to have a realistic theory which provides correct predictions of the behavior 
of a particle similar to quantum theory and, at the same time, respects the description of 
physical reality by EPR as "complete"? Bohm and his followers have attempted a "hidden 
variable theory" , which seemed to satisfy these requirements [6] . The hidden variable theory 
was successfully applied to many different quantum phenomena until 1964, when Bell proved a 
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theorem to show that an inequality, which is violated by certain quantum mechanical statistical 
predictions, can be used to distinguish local hidden variable theory from quantum mechanics 
[7]. Since then, the testing of Bell's inequalities became a standard instrument for the study 
of fundamental problems of quantum theory [8] . The experimental testing of Bell's inequality 
started from the early 1970's. Most of the historical experiments concluded the violation of the 
Bell's inequalities and thus disproved the local hidden variable theory [8] [9] [10] . 

In the following, we examine a simple yet popular realistic model to simulate the behavior of 
the entangled EPR system. This model concerns an ensemble of classically correlated particles 
instead of the quantum mechanical superposition of a particle. In terms of "cats" , this model is 
based on the measurement of a large number of twin cats in which 50% are alive-alive twins and 
50% are dead-dead twins. This model refuses the concept of Schrodinger's cat which requires 
a cat to be alive and dead simultaneously, and each pair of cats involved in a joint detection 
event is in the state of alive-alive and dead-dead simultaneously. 

In this model, we may have three different states: 

(1) State one, each single pair of particles holds defined momenta pi — constant and p 2 = 
constant with pi + p 2 = 0. From pair to pair, the values of p\ and p% may vary significantly. 
The sum of p\ and P2, however, keeps a constant of zero. Thus, each joint detection of the two 
distant particles measures precisely the constant values of p\ and pi and measures p\ + P2 = 0. 
The uncertainties of Api and Ap2 only have statistical meaning in terms of the measurements 
of an ensemble. This model successfully simulated A(p\ + P2) = based on the measurement 
of a large number of classically correlated particle pairs. This is, however, only half of the EPR 
story. Can we have A(xi — X2) = simultaneously in this model? We do have Ax± ~ 00 
and Ax2 ~ 00, otherwise the uncertainty principle will be violated. The position correlation, 
however, can never achieve A(xi — x%) = by any means. 

(2) State two, each single pair of particles holds a well defined position x\ = constant 
and X2 = constant with X\ — X2 — xq. From pair to pair, the values of x\ and X2 may vary 
significantly. The difference of x\ and X2, however, maintains a constant of Xq. Thus, each joint 
detection of the two distant particles measures precisely the constant values of x\ and X2 and 
measures X\ — X2 — Xq. The uncertainties of Aa^ and Ax 2 only have statistical meaning in 
terms of the measurements of an ensemble. This model successfully simulated A(xi — X2) = 
based on the measurement of a large number of classically correlated particle pairs. This is, 
however, only half of the EPR story. Can we have A(pi +P2) = simultaneously in this model? 
We do have Api ~ 00 and Ap2 ~ 00, otherwise the uncertainty principle will be violated. The 
momentum correlation, however, can never achieve A(pi + P2) — by any means. 

The above two models of classically correlated particle pairs can never achieve both A(pi + 
P2) = and A(xi — X2) — 0. What would happen if we combine the two parts together? This 
leads to the third model of classical simulation. 

(3) State three, among a large number of classically correlated particle pairs, we assume 
50% to be in state one and the other 50% state two. The p\ + P2 measurements would have 
50% chance with p\ + p 2 — and 50% chance with p\ + p 2 = random value. On the other 
hand, the X\ — X2 measurements would have 50% chance with x\ — X2 — xq and 50% chance 
with x\ — x 2 = random value. What are the statistical uncertainties on the measurements of 
(pi +P2) and {x\ — X2) in this case? If we focus on only these events of state one, the statistical 
uncertainty on the measurement of (pi +P2) is A(pi +P2) = 0, and if we focus on these events of 
state two, the statistical uncertainty on the measurement of {x\ —X2) is A(x\ —X2) — 0; however, 
if we consider all the measurements together, the statistical uncertainties on the measurements 
of (pi + P2) and (2:1 — X2), are both infinity: A(pi + p%) = 00 and A(xi — X2) = 00. 

In conclusion, classically correlated particle pairs may partially simulate EPR correlation 
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with three types of optimized observations: 



(1) A(pi + p 2 ) - (100%) & A(a;i - x 2 ) = oo (100%); 

(2) A(xi - x 2 ) = (100%) & A(pi + p 2 ) = oo (100%); 

(3) A(pi + pa) = (50%) & A(si - x 2 ) = (50%); 



Within one setup of experimental measurements, only the entangled EPR states result in the 
simultaneous observation of 

A(pi + pa) = (100%) & A(a?i - ar a ) = (100%) 
Api ~ oo, Ap2 ~ oo, Axi ~ oo, Ax 2 ~ oo. 

We thus have a tool, besides the testing of Bell's inequality, to distinguish quantum entangled 
states from classically correlated particle pairs. 

2 Entangled state 

The entangled state of a two-particle system was mathematically formulated by Schrodingcr 
[TT] . Consider a pure state for a system composed of two distinguishable subsystems 

|*) = £>(M)|a>|&> (5) 

a . b 

where {| a)} and {| b)} are two sets of orthogonal vectors for subsystems 1 and 2, respectively. 
If c(a,b) does not factor into a product of the form /(a) x g(b), then it follows that the state 
does not factor into a product state for subsystems 1 and 2: 

p = |*)<*| = £ c(a, b)\a)\b) c*(a', W|(a'| ^ pi x p 2 , (6) 

a, 6 a'.b 7 

where p is the density operator, the state was defined by Schrodinger as an entangled state. 

Following this notation, the first classic entangled state of a two-particle system, the EPR 
state of Eq. {1} and Eq. ([2} , is thus written as: 

\^)epr = ^2 S ( Xl - x z+ x o) \xi)\x 2 ) = ^2 5 (P l + P 2 ) bi )b2 ), (7) 

where we have described the entangled two-particle system as the coherent superposition of the 
momentum eigenstates as well as the coherent superposition of the position eigenstates. The two 
(^-functions in Eq. ([7|) represent, respectively and simultaneously, the perfect position-position 
and momentum-momentum correlation. Although the two distant particles are interaction-free, 
the superposition selects only the eigenstates which are specified by the J-function. We may 
use the following statement to summarize the surprising feature of the EPR state: the values of 
the momentum and the position for neither interaction-free single subsystem is determinated. 
However, if one of the subsystems is measured to be at a certain value of momentum and/or 
position, the momentum and/or position of the other one is 100% determined, despite the 
distance between them. 

It should be emphasized again that Eq. |(7J) is true, simultaneously, in the conjugate space 
of momentum and position. This is different from classically correlated states 

P=^2 5{pi+p 2 )\p 1 )\p 2 ){p 2 \{pi\, (8) 
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or 

P = X! ~ x 2 + xo) \ xi)\x 2 ){x 2 \{xi\. (9) 

Xl,l2 

Eq. ([8]) and Eq. ([9]) represent mixed states. Eq. ([8]) and Eq. ([9]) cannot be true simultaneously as 
we have discussed earlier. Thus, we can distinguish entangled states from classically correlated 
states through the measurements of the EPR inequalities of Eq. (U) . 

Two-photon state of spontaneous parametric down-conversion 

The state of a signal-idler photon pair created in spontaneous parametric down-conversion 
(SPDC) is a typical EPR state [H][I3]- Roughly speaking, the process of SPDC involves 
sending a pump laser beam into a nonlinear material, such as a non-centrosymmetric crystal. 
Occasionally, the nonlinear interaction leads to the annihilation of a high frequency pump 
photon and the simultaneous creation of a pair of lower frequency signal-idler photons forming 
an entangled two-photon state: 

|*) = * 5^*(w, + w 4 -w p )5(k, + ki-k p )a+(k a )ot(k i ) | 0) (10) 

s,i 

where u)j, kj (j = s, i, p) are the frequency and wavevector of the signal (s), idler (i), and pump 
(p), a\ and a\ are creation operators for the signal and the idler photon, respectively, and \&o IS 
the normalization constant. We have assumed a CW monochromatic laser pump, i.e., u p and 
k p are considered as constants. The two delta functions in Eq. (fTO)) are technically named as 
the phase matching condition [12] [14] : 

Lj p = LJ s +L0i, k p = k s + kj. (11) 

The names signal and idler are historical leftovers. The names perhaps came about due to the 
fact that in the early days of SPDC, most of the experiments were done with non-degenerate 
processes. One radiation was in the visible range (and thus easily observable, the signal), while 
the other was in the IR range (usually not measured, the idler). We will see in the following 
discussions that the role of the idler is no any less important than that of the signal. The SPDC 
process is referred to as type-I if the signal and idler photons have identical polarizations, and 
type-II if they have orthogonal polarizations. The process is said to be degenerate if the SPDC 
photon pair has the same free space wavelength (e.g. = X s = 2A p ), and nondegenerate 
otherwise. In general, the pair exit the crystal non-collinearly, that is, propagate to different 
directions defined by the second equation in Eq. (jTTJ) and Snell's law. In addition, the pair may 
also exit collinearly, in the same direction, together with the pump. 

The state of the signal-idler pair can be derived, quantum mechanically, by the first order 
perturbation theory with the help of the nonlinear interaction Hamiltonian. The SPDC inter- 
action arises in a nonlinear crystal driven by a pump laser beam. The polarization, i.e., the 
dipolc moment per unit volume, is given by 

Pi = xf]E j + x^lkEjEk + xfl^EuEt + ... (12) 

where the rath order electrical susceptibility tensor. In SPDC, it is the second or- 

der nonlinear susceptibility x^ that plays the role. The second order nonlinear interaction 
Hamiltonian can be written as 

H = e J^dr xl% EiEjEk (13) 
where the integral is taken over the interaction volume V . 
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It is convenient to use the Fourier representation for the electrical fields in Eq. (JTJ 

E(r, t) = J dk[ E (-)(k)e- l(w(k)t - k - r) + E(+)(k)e^ k )*- k - r > ]. (14) 

Substituting Eq. (fT4")) into Eq. (flU)) and keeping only the terms of interest, we obtain the SPDC 
Hamiltonian in the interaction representation: 

H lnt {t) (15) 



J dv J dk.dk, ^E^e'^^ 



where h.c. stands for Hermitian conjugate. To simplify the calculation, we have also assumed 
the pump field to be a monochromatic plane wave with wave vector k p and frequency u p . 

It is easily noticeable that in Eq. (fT5"|) , the volume integration can be done for some simplified 
cases. At this point, we assume that V is infinitely large. Later, we will see that the finite size 
of V in longitudinal and/or transversal directions may have to be taken into account. For an 
infinite volume V , the interaction Hamiltonian Eq. (|15p is written as 



H mt {t) =e J dk s dk, xZl EffE^EtJ 6(k p k s i^K-^)-^))* + h . c . (i 6 ) 

It is reasonable to consider the pump field to be classical, which is usually a laser beam, and 
quantize the signal and idler fields, which are both at the single-photon level: 



£(-)(k) = i^at(k), £M(k) = i^afk), (17) 

where a^(k) and a(k) are photon creation and annihilation operators, respectively. The state 
of the emitted photon pair can be calculated by applying the first order perturbation 

m = -^JdtH mt (t)\0). (18) 

By using vacuum |0) for the initial state in Eq. (I18|) . we assume that there is no input radiation 
in any signal and idler modes, that is, we have a spontaneous parametric down conversion 
(SPDC) process. 

Further assuming an infinite interaction time, evaluating the time integral in Eq. (| 18|> and 
omitting altogether the constants and slow (square root) functions of u>, we obtain the entangled 
two-photon state of Eq. (jTUJ) in the form of an integral [T3] : 

|*) = * /<*A% - to s {k s ) - Wi(ki)]5(k p - k s - k0at(k,)aj(ki)|0) (19) 

where is a normalization constant which has absorbed all omitted constants. 



The way of achieving phase matching, i.e., the delta functions, in Eq. (fl"9]l basically deter- 
mines how the signal-idler pair "looks" . For example, in a negative uniaxial crystal, one can 
use a linearly polarized pump laser beam as an extraordinary ray of the crystal to generate a 
signal-idler pair both polarized as the ordinary rays of the crystal, which is defined as type-I 
phase matching. One can alternatively generate a signal-idler pair with one ordinary polarized 
and another extraordinary polarized, which is defined as type II phase matching. Fig. [1] shows 
three examples of an SPDC two-photon source. All three schemes have been widely used for 
different experimental purposes. Technical details can be found in text books and research 
references in nonlinear optics. 

The two-photon state in the forms of Eq. (fT0|) or Eq. (fT9|) is a pure state, which math- 
ematically describes the behavior of a signal-idler photon pair. The surprise comes from the 
coherent superposition of the two-photon modes: 
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Figure 1: Three widely used SPDC setups, (a) Type-I SPDC. (b) Collinear degenerate 
type-II SPDC. Two rings overlap at one region, (c) Non-collinear degenerate type-II SPDC. 
For clarity, only two degenerate rings, one for e-polarization and the other for o-polarization, 
are shown. Notice, the color rainbows represent the distribution function of a signal-idler 
pair. One signal-idler pair yields the entire rainbow. 



Does the signal or the idler photon in the EPR state of Eq. (jTUJ) or Eq. (fl9)) have a defined 
energy and momentum regardless of whether we measure it or not? Quantum mechanics 
answers: No! However, if one of the subsystems is measured with a certain energy and 
momentum, the other one is determined with certainty, despite the distance between them. 



It is indeed a mystery from a classical point of view. There has been, nevertheless, classical 
models to avoid the surprises. One of the classical realistic models insists that the state of 
Eq. (JXOJ) or Eq. (|19p only describes the behavior of an ensemble of photon pairs. In this model, 
the energy and momentum of the signal photon and the idler photon in each individual pair are 
defined with certain values and the resulting state is a statistical mixture. Mathematically, it 
is incorrect to use a pure state to characterize a statistical mixture. The concerned statistical 
ensemble should be characterized by the following density operator 

dk s dkiS(u! p - u; s - uJi)6(k p - k s - k») «4(k s ) a|(kj) | 0)(0 | a s (k s ) aj(ki) (20) 

which is very different from the pure state of SPDC. We will show later that a statistical mixture 
of Eq. (|20[) can never have delta- function- like two-photon temporal and/or spatial correlation 
that is shown by the measurement of SPDC. 

For finite dimensions of the nonlinear interaction region, the entangled two-photon state of 
SPDC may have to be estimated in a more general format. Following the earlier discussions, 
we write the state of the signal-idler photon pair as 



/dk«dkiF(k„k i )aJ(k,)ot(k i )|Q) (21) 



where 



F(k a ,kj) = eS(uj p - lj s - LUt) f(A z L) h tr {Ki + k 2 ) 
f(A z L) = ( dze- i{k "- k "- k " )z 



L 



M*i + *2)= / dphUpje-^'+^-P (22) 



A 

A z = kp k s 



where e is named as the parametric gain index, e is proportional to the second order electric sus- 
ceptibility an d is usually treated as a constant, L is the length of the nonlinear interaction, 
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the integral in k is evaluated over the cross section A of the nonlinear material illuminated by 
the pump, p is the transverse coordinate vector, itj (with j = s, i) is the transverse wavevector 
of the signal and idler, and /(| p |) is the transverse profile of the pump, which can be treated 
as a Gaussion in most of the experimental conditions. The functions f(A z L) and ht r (Ki + K2) 
turn to (5-functions for an infinitely long (L ~ 00) and wide (A ~ 00) nonlinear interaction 
region. The reason we have chosen the form of Eq. (|22|) is to separate the "longitudinal" and 
the "transverse" correlations. We will show that 5{lo v — uo s — u>i) and f(A z L) together can be 
rewritten as a function of u> s — u>i. To simplify the mathematics, we assume near co-linearly 
SPDC. In this situation, | n s ,i | <C | k Sj i |. 



Basically, the function f(A z L) determines the "longitudinal" space-time correlation. Find- 
ing the solution of the integral is straightforward: 



f{A z L) = [ dze- l{k "- k "- k ">- = e-^ L / 2 sinc{A z L/2). 
Jo 



(23) 



Now, we consider f(A z L) with 6(lu p — uj s — cjj) together, and taking advantage of the 
(5-function in frequencies by introducing a detuning frequency VL to evaluate function f(A z L): 

Ui = lu°-Q (24) 
Q = (lj s - L0i)/2. 

The dispersion relation k{oj) allows us to express the wave numbers through the frequency 
detuning fl: 



dk 



k s « k(io° s ) + n— =k(u° s ) + -, 



duj 
dk 



n 



h a fc( w °)-fi^ = fc( w °)_- (25) 
dio ui° Ui 

where u s and ui are group velocities for the signal and the idler, respectively. Now, we connect 
A z with the detuning frequency f2: 



A z — k n k 



sz <H2 



k p ~ ^(k s f - (K s f - ^(hf - (K^ 
(Zs) 2 | (ft) 2 

2 a 2 ki 



~ k P -k s -k i + ^ + ^- (26) 



k - fc(^) - k(u>») + 

U s U 



. n n , (k s ) 2 , (k 4 ) 2 



2/bc 2/bl 



where D = l/u s — 1/ui. We have also applied k p — k(ui®) — k(cu®) — and | K St i | <C | k s 
The "longitudinal" wavevector correlation function is rewritten as a function of the detuning 
frequency = (uj s — o->i)/2: f(A z L) = f(flDL). In addition to the above approximations, 
we have inexplicitly assumed the angular independence of the wavevector k = n(9)uj/c. For 
type II SPDC, the refraction index of the extraordinary-ray depends on the angle between the 
wavevector and the optical axis and an additional term appears in the expansion. Making the 
approximation valid, we have restricted our calculation to a near-collinear process. Thus, for a 
good approximation, in the near-collinear experimental setup 

A Z L = flDL = (u s - u)i)DL/2. (27) 
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Type-I degenerate SPDC is a special case. Due to the fact that u s — Uj, and hence, D = 0, 
the expansion of k(u>) should be carried out up to the second order. Instead of (|2"T|) . we have 



A Z L ^ -Q 2 D'L = -{uj s - uj t ) 2 D'L/4: 



(28) 



where 




,0 



The two-photon state of the signal-idler pair is then approximated as 




dfl dn s (Mi /(Q) h tr (n s + K>i) al(uj® + f2, k s ) a\ (id® — ft, Ki)\ ) 



(29) 



where the normalization constant has been absorbed into /(^)- 

3 Correlation measurement of entangled state 

EPR state is a pure state which characterizes the behavior of a pair of entangled particles. 
In principle, one EPR pair contains all information of the correlation. A question naturally 
arises: Can we then observe the EPR correlation from the measurement of one EPR pair? 
The answer is no. Generally speaking, we may never learn any meaningful physics from the 
measurement of one particle or one pair of particles. To learn the correlation, an ensemble 
measurement of a large number of identical pairs are necessary, where "identical" means that 
all pairs which are involved in the ensemble measurement must be prepared in the same state, 
except for an overall phase factor. This is a basic requirement of quantum measurement theory. 

Correlation measurements are typically statistical and involve a large number of measure- 
ments of individual quanta. Quantum mechanics does not predict a precise outcome for a 
measurement. Rather, quantum mechanics predicts the probabilities for certain outcomes. In 
photon counting measurements, the outcome of a measurement is either a "yes" (a count or a 
"click") or a "no" (no count). In a joint measurement of two photon counting detectors, the 
outcome of "yes" means a "yes-yes" or a "click-click" joint registration. If the outcome of a 
joint measurement shows 100% "yes" for a certain set of values of a physical observable or a 
certain relationship between physical variables, the measured quantum system is correlated in 
that observable. As a good example, EPR's gedankenexperiment suggested to us a system of 
quanta with perfect correlation S(x\ — X2 + £o) in position. To examine the EPR correlation, we 
need to have a 100% "yes" when the positions of the two distant detectors satisfy x\ — X2 = xo, 
and 100% "no" otherwise, when x\ — X2 f= xq. To show this experimentally, a realistic approach 
is to measure the correlation function of \ f{x\ — X2)\ 2 by observing the joint detection counting 
rates of Ri t 2 oc \f(x\ — X2)\ 2 while scanning all possible values of x\ — x^. In quantum optics, 
this means the measurement of the second-order correlation function, or (ri, tx\ ?2, ^2), in 
the form of longitudinal correlation G^{t\ — T2) and/or transverse correlation G^(pi — P2), 
where Tj = tj — Zj/c, j — 1,2, and pj is the transverse coordinate of the jth point-like photon 
counting detector. 

Now, we study the two-photon correlation of the entangled photon pair of SPDC. The 
probability of jointly detecting the signal and idler at space-time points (ri,ii) and (^,£2) is 
given by the Glauber theory [15] : 



the detection events at space-time points (ri,ii) and (r2,t2). The expectation value of the joint 
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detection operator is calculated by averaging over the quantum states of the signal-idler photon 
pair. For the two-photon state of SPDC, 

G«(r x , t i; r a , t 2 ) = | (0 | £(+>(r a , t 2 )£ (+) (ri, h) | * ) | 2 = | ^(r X) i x ; r 2 , f a ) | 2 (31) 

where ^ ) is the two-photon state, and \I/(ri , t\\ r 2 , i 2 ) is named the effective two-photon wave- 
function. To evaluate G^(ri, t\; r 2 , t%) and ip(ri 7 1\; r 2 , t 2 ), we need to propagate the field 
operators from the two-photon source to space-time points (ri,ti) and (r a ,i a ). 

In general, the field operator E^(r,t) at space-time point (r, t) can be written in terms 
of the Green's function, which propagates a quantized mode from space-time point (ro,£o) to 

(r,t) mm- 

£«(r,i) = £ ff(k,r-r ,i-t )-B (+) (k,ro,to). (32) 

k 

where g(k, r— ro, t— to) is the Green's function, which is also named the optical transfer function. 
For a different experimental setup, g(k, r — ro,t — to) can be quite different. To simplify the 
notation, we have assumed one polarization. 

Di 




D 2 



Figure 2: Collinear propagated signal-idler photon pair, either degenerate or non-degenerate, 
are received by two distant point photo-detectors D\ and D 2 , respectively, for longitudinal 
(n — t 2 ) and transverse G^ 2 ' (pi — p 2 ) correlation measurements. To simplify the math- 
ematics, we assume paraxial approximation is applicable to the signal-idler fields. The z\ 
and z 2 are chosen along the central wavevector k° and k° . 

Considering an idealized simple experimental setup, shown in Fig. [2l in which collinear 
propagated signal and idler pairs are received by two point photon counting detectors D\ and 
Z) 2 , respectively, for longitudinal G^ 2 \t\ — t 2 ) and transverse G^ 2 \pi — p 2 ) correlation measure- 
ments. To simplify the mathematics, we further assume paraxial experimental condition. It is 
convenient, in the discussion of longitudinal and transverse correlation measurements, to write 
the field (rj, tj) in terms of its longitudinal and transversal space-time variables under the 
Fresnel paraxial approximation: 

£(+)(/>,, z^) (33) 

— J dui dti g(K,,Lu; pj, zj) e~ luJtj a(u;, k) = J du> dn^f(K,uj; pj, Zj) e~ luJTj a(uj, k) 

where g(n, u; pj, Zj) — 7(k, pj, Zj)e luJZj / c is the spatial part of the Green's function, pj and 
Zj are the transverse and longitudinal coordinates of the jth photo-detector and k is the trans- 
verse wavevector. We have chosen z = and t = at the output plane of the SPDC. For 
convenience, all constants associated with the field are absorbed into <?(k, lu] pj, Zj). 
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The two-photon effective wavefunction ^(pi, z%, t%; p 2 , Z2, £2) is thus calculated as follows 



&(p 1 ,Zl,ti;p2,Z2,t 2 ) 

= (0 / du' dn' <?(>?', W; p2, z 2 ) e ~ IUJ * 2 a(w, k') 



x J duj" dti" g(n" : a>"; pi, zi) e la, " t:L a(a;", ft') 

x y c?r2 d/? s dKif(fl) h tr (K s + Ki) al(iu® + k s ) a\(uj® ~ Cl 7 Ki)\0) 

x J dQ dit s dKi /(O) /^(ft, + ft) e-^( T ^ T2 ) 7 (K s , O; fa, Zl ) 7 (ft, -fi; p 2 , z 2 ). (34) 

Although Eq. (j3"4")) cannot be factorized into a trivial product of longitudinal and transverse 
integrals, it is not difficult to measure the temporal correlation and the transverse correlation 
separately by choosing suitable experimental conditions. 

Experiments may be designed for measuring either temporal (longitudinal) or spatial (trans- 
verse) correlation only. Thus, based on different experimental setups, we may simplify the 
calculation to either the temporal (longitudinal) part: 

*(Ti;T 2 ) = for ,(u ° Tl+u - T2) J dO/(0)e- ,!!(Tl - T2) = * e" l(u ° Tl+ ' J?T2, f Tl - T2 {/(l])} (35) 

or the spatial part: 

^(p 1 ,z 1 ;p 2 ,z 2 ) = *o / dK s dKih tr (K s + ft) g(K s ,uj s ; pi, z x ) s(ft,Wj; p%,zi). (36) 



In Eq. (f35|) . J- Tl - T2 {/(£!)} is the Fourier transform of the spectrum amplitude function /(f2). 
In Eq. (|36"1) , we may treat /i{ r (fti + ft) ~ <5(ft + ft) by assuming certain experimental conditions. 



Two-photon temporal correlation 

To measure the two-photon temporal correlation of SPDC, we select a pair of transverse 
wavevectors hi a — in Eq. (|34|) by using appropriate optical apertures. The effective two- 
photon wavefunction is thus simplified to that of Eq. (I35|) 

*(ri;r 2 ) S * e- ,; ^ Tl+ ^ T2 ) J dttf{VL)e- rn{Tl - T ^ (37) 

where, again, J- T1 _ T2 j/(f2)} is the Fourier transform of the spectrum amplitude function /(il). 
Eq. (|37[) indicates a 2-D wavepacket: a narrow envelope along the T\ — r 2 axis with constant 
amplitude along the T\ + r 2 axis. In certain experimental conditions, the function /(f2) of 
SPDC can be treated as constant from — 00 to 00 and thus T Ti -t 2 ~ <5( T i — T 2)- I n this case, 
for fixed positions of D\ and D 2 , the 2-D wavepacket means the following: the signal- idler 
pair may be jointly detected at any time; however, if the signal is registered at a certain time 
t\, the idler must be registered at a unique time of t 2 ~ t\ — (zi — z 2 )/c. In other words, 
although the joint detection of the pair may happen at any times of t\ and t 2 with equal 
probability (A(ti -I- t 2 ) ~ 00), the registration time difference of the pair must be a constant 
A(ti — i 2 ) ~ 0. A schematic of the two-photon wavepacket is shown in Fig. [3l It is a non- 
factorizeable 2-D wavefunction indicating the entangled nature of the two-photon state. The 
longitudinal correlation function G' 2 '(ti — t 2 ) is thus 

G^i-r 2 )oc|^_ T2 {/(ft)}| 2 , 
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Figure 3: A schematic envelope of a two-photon wavepacket with a Gaussian shape along 
n — T2 corresponding to a Gaussian function of /(ft). In the case of SPDC, the envelope 
is close to a (5-function in t% — t-i corresponding to a broad-band /(ft) = constant. The 
wavepacket is uniformly distributed along n + T2 due to the assumption of lu p — constant. 

which is a ^-function-like function in the case of SPDC. Thus, we have shown the entangled 
signal-idler photon pair of SPDC hold a typical EPR correlation in energy and time: 



Now we examine a statistical model of SPDC for temporal correlation. As we have discussed 
earlier, realistic statistical models have been proposed to simulate the EPR two-particle state. 
Recall that for a mixed state in the form of 



where Pj is the probability for specifying a given set of state vectors the second-order 

correlation function of fields E(ri,ti) and £^2^2) is given by 



A(w s + uJi) - & A(ti - h) ~ 
with Aw s ~ 00, Auji ~ 00, Afi ~ 00, Ai 2 ~ 00. 



p=E-pji*i><*ii 



= Tr[^^-J(ri,ti) £(-)(r 2 ,t 2 ) £«(r 2 ,t 2 ) £(+>(n,ti) ] 

= ^ Pj ( % I £(-)( ri , ti) (r 2 , t 2 ) # (+) (r 2 , *a) # (+) fa, <i) | * rf ) 



which is a weighted sum over all individual contributions of Gj . Considering the following 
simplified version of Eq. (|20p to simulate the state of SPDC as a mixed state: 




dft I /(ft) | 2 a^{u° s + ft) a\u)l - ft)| 0)<0 |a(o;? - ft) a(w s ° + ft), (38) 



with 



l*n> 



a 



t (w s ° + ft)a t (w 4 -ft)|0), =dft|/(ft)| 2 . 



(39) 
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It is easy to find Gq (t\ — t 2 ) = constant, and thus (ti—t 2 ) = constant. This means that the 
uncertainty of the measurement on t\—t 2 for the mixed state of Eq. (|38|) is infinite: A(i\ — t 2 ) ~ 

00. Although the energy (frequency) or momentum (wavevector) for each photon may be defined 
with constant values pair by pair, the corresponding temporal correlation measurement of the 
ensemble can never achieve a <5-function-like relationship. In fact, the correlation is undefined, 

1. e., taking an infinite uncertainty. Thus, the statistical model of SPDC cannot satisfy the EPR 
inequalities of Eq. (|4]) . 

Two-photon spatial correlation 

Similar to that of the two-photon temporal correlation, as an example, we analyze the 
effective two-photon wavefunction of the signal-idler pair of SPDC. To emphasize the spatial 
part of the two-photon correlation, we choose a pair of frequencies w a and loi with lu s +u>i — u> p . 
In this case, the effective two-photon wavefunction of Eq. ([34| is simplified to that of Eq. ([36]) 

^(pi,zi; p 2 ,z 2 ) = * J dn s dR i S(U s + K i )g(K s ,uJ s ,pi,Zi)g(Ki,uj i ,p 2 ,Z2) 

where we have assumed h tr (H s + Hi) ~ 5(H S + Hi), which is reasonable by assuming a large 
enough transverse cross-session laser beam of pump. 

We now design a simple joint detection measurement between two point photon counting 
detectors D\ and D 2 located at (pi, zi) and (p 2 , z 2 ), respectively, for the detection of the signal 
and idler photons. We have assumed that the two-photon source has a finite but large transverse 
dimension. Under this simple experimental setup, the Green's function, or the optical transfer 
function describing arm-j, j = 1, 2, in which the signal and the idler freely propagate to photo- 
detector Di and D 2 , respectively, is given by Eq. jA— 5D of the Appendix. Substitute the 
9j(oJ, k;zj, pj), j — 1,2, into Eq. (|36p . the effective wavefunction is then given by 

*(pi,zi;p 2 ,z 2 ) (40) 

x / dpsdpiGQpi-psl^e^GUfe-pil—)^^ 
Ja czi cz 2 

where p s (k s ) and pi (Hi) are the transverse coordinates (wavevectors) for the signal and the 
idler fields, respectively, defined on the output plane of the two-photon source. The integral of 
dp s and dpi is over area A, which is determined by the transverse dimension of the nonlinear 
interaction. The Gaussian function G(\a\,(5) = e 1 ^/ 2 ""' represents the Fresnel phase factor 
that is defined in the Appendix. The integral of dn s and di?i can be evaluated easily with the 
help of the EPR type two-phonon transverse wavevector distribution function 6(k s + «i): 

J dn s dKi 5(k s + Hi) e iR ^ e^* ~ 6(p s - p % ). (41) 

Thus, we have shown that the entangled signal-idler photon pair of SPDC holds a typical 
EPR correlation in transverse momentum and position while the correlation measurement is on 
the output plane of the two-photon source, which is very close to the original proposal of EPR: 

A(H S + Hi) - & A(p s - p^ ~ 
with Ak s ~ 00, A.Hi ~ 00, Ap s ~ 00, Api ~ 00. 

In EPR's language, we may never know where the signal photon and the idler photon are 
emitted from the output plane of the source. However, if the signal (idler) is found at a certain 
position, the idler (signal) must be observed at a corresponding unique position. The signal and 
the idler may have also any transverse momentum. However, if the transverse momentum of the 
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signal (idler) is measured at a certain value in a certain direction, the idler (signal) must be of 
equal value but pointed to a certain opposite direction. In collinear SPDC, the signal-idler pair 
is always emitted from the same point in the output plane of the two-photon source, p s = pi, 
and if one of them propagates slightly off from the collinear axes, the other one must propagate 
to the opposite direction with 

The interaction of spontaneous parametric down-conversion is nevertheless a local phe- 
nomenon. The nonlinear interaction coherently creates mode-pairs that satisfy the phase match- 
ing conditions of Eq. which are also named as energy and momentum conservation. The 
signal-idler photon pair can be excited to any of these coupled modes or in all of these coupled 
modes simultaneously, resulting in a particular two-photon superposition. It is this superposi- 
tion among those particular "selected" two-photon states which allows the signal-idler pair to 
come out from the same point of the source and propagate to opposite directions with k s = 

The two-photon superposition becomes more interesting when the signal-idler is separated 
and propagated to a large distance, either by free propagation or guided by optical components 
such as a lens. A classical picture would consider the signal photon and the idler photon 
independent whenever the pair is released from the two-photon source because there is no 
interaction between the distant photons in free space. Therefore, the signal photon and the 
idler photon should have independent and random distributions in terms of their transverse 
position pi and p2- This classical picture, however, is incorrect. It is found that the signal- 
idler two-photon system would not lose its entangled nature in the transverse position. This 
interesting behavior has been experimentally observed in quantum imaging by means of an EPR 
type correlation in transverse position. The sub-diffraction limit spatial resolution observed in 
the "quantum lithography" experiment and the nonlocal correlation observed in the "ghost 
imaging" experiment are both the results of this peculiar superposition among those "selected" 
two-photon amplitudes, namely that of two-photon superposition, corresponding to different 
yet indistinguishable alternative ways of triggering a joint photo-electron event at a distance. 
Two-photon superposition does occur in a distant joint detection event of a signal-idler photon 
pair. There is no surprise that one has difficulties facing this phenomenon. The two-photon 
superposition is a nonlocal concept in this case. There is no counterpart for such a concept in 
classical theory and it may never be understood classically. 

Now we consider propagating the signal- idler pair away from the source to (p\,Z\) and 
{P2,Z2)i respectively, and taking the result of Eq. (|4ip . i.e., p s = Pi = Po on the output plane 
of the SPDC source, the effective two-photon wavefunction becomes 

*(pi,zi;p 2 ,z 2 ) (42) 



(2TTC) 2 ZlZ2 J a cz 1 cz 2 

where po is defined on the output plane of the two-photon source. Eq. ([%2l indicates that the 
propagation-diffraction of the signal and the idler cannot be considered as independent. The 
signal-idler photon pair are created and diffracted together in a peculiar entangled manner. 
This point turns out to be both interesting and useful when the two photodetectors coincided, 
or are replaced by a two-photon sensitive material. Taking z\ = z-i and p\ = p 2 , Eq. (|42[) 
becomes 

*(p,z;p,z) = --^^e^ J^ d p G(\p-p \,^) (43) 

where uj p is the pump frequency, which means that the signal-idler pair is diffracted as if 
they have twice the frequency or half the wavelength. This effect is named as "two-photon 
diffraction" . This effect is useful for enhancing the spatial resolution of imaging. 
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4 Quantum imaging 



Although questions regarding fundamental issues of quantum theory still exist, quantum 
entanglement has started to play important roles in practical engineering applications. Quan- 
tum imaging is one of these exciting areas [18) . Taking advantage of entangled states, Quantum 
imaging has so far demonstrated two peculiar features: (1) enhancing the spatial resolution of 
imaging beyond the diffraction limit, and (2) reproducing ghost images in a "nonlocal" man- 
ner. Both the apparent "violation" of the uncertainty principle and the "nonlocal" behavior 
of the momentnm-momentum position-position correlation are due to the two-photon coherent 
effect of entangled states, which involves the superposition of two-photon amplitudes, a non- 
classical entity corresponding to different yet indistinguishable alternative ways of triggering a 
joint-detection event in the quantum theory of photodetection. In this section, we will focus 
our discussion on the physics of imaging resolution enhancement. The nonlocal phenomenon of 
ghost imaging will be discussed in the following section. 

The concept of imaging is well defined in classical optics. Fig. [4] schematically illustrates a 
standard imaging setup. A lens of finite size is used to image the object onto an image plane 




Figure 4: A lens produces an image of an object in the plane defined by the Gaussian thin 
lens equation l/si + l/s = 1//. The concept of an image is based on the existence of a 
point-to-point relationship between the object plane and the image plane. 



which is defined by the "Gaussian thin lens equation" 

f + f-i (44) 

Si S J 

where s is the distance between object and lens, / is the focal length of the lens, and Sj is the 
distance between the lens and image plane. If light always follows the laws of geometrical optics, 
the image plane and the object plane would have a perfect point-to-point correspondence, which 
means a perfect image of the object, either magnified or demagnified. Mathematically, a perfect 
image is the result of a convolution of the object distribution function f(p ) and a (5-function. 
The (5-function characterizes the perfect point-to-point relationship between the object plane 
and the image plane: 

F(pi)= [ dp f(p )5(p + ^) = .f(p )®5(p + ^) (45) 
Job] m rn 

where p and pi are 2-D vectors of the transverse coordinate in the object plane and the image 
plane, respectively, and to is the magnification factor. The symbol ® means convolution. 

Unfortunately, light behaves like a wave. The diffraction effect turns the point-to-point 
correspondence into a point-to- "spot" relationship. The 5-function in the convolution of Eq. (|45|) 
will be replaced by a point-spread function. 

F{p i )= f dp f(p ) S omb[--\p + ^\]=f(p )® S omb[--\p + ^\] (46) 

Job] S C' TO S C 1 TO 
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where 

sombyx) = , 

x 

and J\(x) is the first-order Bessel function, R is the radius of the imaging lens. R/s a is named as 
the numerical aperture of the imaging system. The finite size of the spot, which is defined by the 
point-spread function, determines the spatial resolution of the imaging setup, and thus, limits 
the ability of making demagnified images. It is clear from Eq. (|46p . the use of a larger imaging 
lens and shorter wavelength light of source will result in a narrower point-spead function. To 
improve the spatial resolution, one of the efforts in the lithography industry is the use of shorter 
wavelengths. This effort is, however, limited to a certain level because of the inability of lenses 
to effectively work beyond a certain "cutoff" wavelength. 




Figure 5: Typical imaging setup. A lens of finite size is used to produce a demagnified image 
of a object with limited spatial resolution. Replacing classical light with an entangled N- 
photon system, the spatial resolution can be improved by a factor of N, despite the Rayleigh 
diffraction limit. 

Eq. (|46[) imposes a diffraction limited spatial resolution on an imaging system while the 
aperture size of the imaging system and the wavelength of the light source are both fixed. This 
limit is fundamental in both classical optics and in quantum mechanics. Any violation would 
be considered as a violation of the uncertainty principle. 

Surprisingly, the use of quantum entangled states gives a different result: by replacing 
classical light sources in Fig. [5] with entangled N-photon states, the spatial resolution of the 
image can be improved by a factor of N, despite the Rayleigh diffraction limit. Is this a violation 
of the uncertainty principle? The answer is no! The uncertainty relation for an entangled N- 
particle system is radically different from that of N independent particles. In terms of the 
terminology of imaging, what we have found is that the somb(x) in the convolution of Eq. (|46|) 
has a different form in the case of an entangled state. For example, an entangled two-photon 
system has 

R 2u> , ^ pi , 
x = \Po H • 

s c m 

Comparing with Eq. (|46|) . the factor of 2oj yields a point-spread function half the width of that 
from Eq. (|46|) and results in a doubling spatial resolution for imaging. 

It should be further emphasized that one must not confuse a "projection" with an image. 
A projection is the shadow of an object, which is obviously different from the image of an 
object. Fig. [6] distinguishes a projection shadow from an image. In a projection, the object- 
shadow correspondence is essentially a "momentum" correspondence, which is defined only by 
the propagation direction of the light rays. 



We now analyze classical imaging. The analysis starts with the propagation of the field 
from the object plane to the image plane. In classical optics, such propagation is described by 
an optical transfer function h(r — r , t — t ), which accounts for the propagation of all modes of 
the field. To be consistent with quantum optics calculations, we prefer to work with the single- 
mode propagator <?(k, r — r , i — £q), and to write the field E(r,t) in terms of its longitudinal 
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Figure 6: Projection: a light source illuminates an object and no image forming system is 
present, no image plane is denned, and only projections, or shadows, of the object can be 
observed. 



(z) and transverse (p) coordinates under the Fresnel paraxial approximation: 

E(p,z,t)= \ dudn £(re, w) 5 (re, w; p, z) e -1 "* (47) 



where E(oj, re) is the complex amplitude of frequency to and transverse wave- vector re. In Eq. (I47[) 
we have taken zq = and to = at the object plane as usual. To simplify the notation, we 
have assumed one polarization. 

Based on the experimental setup of Fig. [5j 17 (re, u; p, z) is found to be 



g(K,w;pi, s + Si) 

dp f dp l {A(p )e™-?°}{ 



obj 



lens 



o G V\Pl ~P< . 

lire s Q cs 



CS„ i 



2cs. 



)} 



(48) 



where p a , pi, and f>i are two-dimensional vectors defined, respectively, on the object, the lens, 
and the image planes. The first curly bracket includes the object-aperture function A{p Q ) and 
the phase factor e %K contributed to the object plane by each transverse mode re. Here we 
have assumed a far-field finite size source. Thus, a phase factor e tK ' Po appears on the object 
plane of z = 0. If a collimated laser beam is used, this phase factor turns out to be a constant. 
The terms in the second and the fourth curly brackets describe free-space Fresnel propagation- 
diffraction from the source/object plane to the imaging lens, and from the imaging lens to 
the detection plane, respectively. The Fresnel propagator includes a spherical wave function 
e^fe-^VOj ~ z k) and a Fresnel phase factor G(\a\,0) = e 1 ^/ 2 ^ 2 = eMA-fl.lVM**-**). 
The third curly bracket adds the phase factor, G(\pi\, — ^4) = e~ % ^r ; which is introduced by 
the imaging lens. 

Applying the properties of the Gaussian function, Eq. (j48|) can be simplified into the fol- 
lowing form 



g(re, uj;pi,z = s + s») 
, ,2 



{2-KC) 2 SoSi 



a if (s +s 



cm,—) 



dp A{p )G{\p l—)e iK ?° 



obj 



CS C 



x / dpiG(\pi\ 

' lens 



2,± 

C So 



Si J 



(49) 



The image plane is defined by the Gaussian thin- lens equation of Eq. (|4"4")l . Hence, the second 
integral in Eq. (f4"9"l) simplifies and gives, for a finite sized lens of radius R, the so called point- 
spread function of the imaging system: somb(x) = 2J\{x)/x, where x = [— —\p a + Pi/m\], J\{x) 
is the first-order Bessel function and m = Si/s Q is the magnification of the imaging system. 
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Substituting the result of Eqs. ([19"]) into Eq. (|TF]) enables one to obtain the classical self- 
correlation of the field, or, cquivalently, the intensity on the image plane 



I(pi,Zi,ti) = (E*(p i ,z i ,ti)E(p i ,z i ,ti)) (50) 

where (...) denotes an ensemble average. We assume monochromatic light for classical imaging 
as usual. |3 

Case (I): incoherent imaging. The ensemble average of (E*(k,lu) E(k' ,oj)) yields zeros 
except when ft = k' . The image is thus 

I(pi) oc [ dp \A(p )\ 2 \somb{——\p + —\]\ 2 . (51) 



An incoherent image, magnified by a factor of to, is thus given by the convolution between 
the squared moduli of the object aperture function and the point-spread function. The spatial 
resolution of the image is thus determined by the finite width of the |som6| 2 -function. 

Case (II): coherent imaging. The coherent superposition of the k modes in both E*(pi,r) 
and E(pi,r) results in a wavepacket. The image, or the intensity distribution on the image 
plane, is thus 

/(#)« I dp A(p )e t ^^ 2 somb[--\p + ^\} 2 . (52) 
J obj s Q c to 

A coherent image, magnified by a factor of to, is thus given by the squared modulus of the 
convolution between the object aperture function (multiplied by a Fresnel phase factor) and 
the point-spread function. 

For Si < s and s > /, both Eqs. (fSTj) and ([52]) describe a real demagnified inverted image. 
In both cases, a narrower soTO&-function yields a higher spatial resolution. Thus, the use of 
shorter wavelengths allows for improvement of the spatial resolution of an imaging system. 



To demonstrate the working principle of quantum imaging, we replace classical light with 
an entangled two-photon source such as spontaneous parametric down-conversion (SPDC) and 
replace the ordinary film with a two-photon absorber, which is sensitive to two-photon transition 
only, on the image plane. We will show that, in the same experimental setup of Fig. [5j an 
entangled two-photon system gives rise, on a two-photon absorber, to a point-spread function 
half the width of the one obtained in classical imaging at the same wavelength. Then, without 
employing shorter wavelengths, entangled two-photon states improve the spatial resolution of 
a two-photon image by a factor of 2 [19] [20] . We will also show that the entangled two-photon 
system yields a peculiar Fourier transform function as if it is produced by a light source with 
A/2. 

In order to cover two different measurements, one on the image plane and one on the 
Fourier transform plane, we generalize the Green's function of Eq. (pf8|) from the image plane 
of z — s Q + Si to an arbitrary plane of z = s a + d, where d may take any values for different 
experimental setups: 

g{Kj,ujj;p k .,z = s + d) 

dp Q I dpi A(p ) {-^- e^-Po e i^fs _ ^ 

obj J lens ZTTCS CS d 

xG(|p ; |,-^){-^ e ^G(|p1-p fc |,^)}, (53) 
cj Zircd cd 



2 Even if assuming a perfect lens without chromatic aberration, Fresnel diffraction is wavelength dependent. 
Hence, large broadband (Au> ~ oo) would result in blurred images in classical imaging. Surprisingly, the situation 
is different in quantum imaging: no aberration blurring. 
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where p , pi, and pj are two-dimensional vectors defined, respectively, on the (transverse) output 
plane of the source (which coincide with the object plane), on the transverse plane of the imaging 
lens and on the detection plane; and j ~ s,i, labels the signal and the idler; k = 1,2, labels 
the photodetector D\ and The function A(p ) is the object-aperture function, while the 
terms in the first and second curly brackets of Eq. (|53|) describe, respectively, free propagation 
from the output plane of the source/object to the imaging lens, and from the imaging lens to 
the detection plane. 

Similar to the earlier calculation, by employing the second and third expressions given in 
Eq. ( pT^l ), Eq. ([53]) simplifies to 

g{Kj,Uj\p k ,z = s + d) 

= jf^ , ^ {s ° +d) G{\h\^) I dp A(p )G(\p \^)e^° 
(2iTcYs d cd J obl cs 

xf dp lG m^A-\])e-^^. (54) 

Jlens C S CI J 



Substituting the Green's functions into Eq. (|34|) . the effective two-photon wavefunction 
\&(pi, z; p2, z) is thus 

y(p u z-p 2 ,z) = * fdSlf(n)G(\^\,^)GQp2\A) 
/ cd cd 



dp A(p ) GQp a \, -L) / dp' a A(p' ) G(\p' \, -L) 

obj CS Q J Q ijj CS 

x / dpiGm^ + l-^e-^^y* 

Jlens C S a J 

x / d? l G{\?M-l- + i-k)e-** l » + * >Pt 

Jlens C S CI J 

x / dn s dKi S(k s + Hi) e^.-A+fc-p'.) (55) 



where we have absorbed all constants into "fo, including the phase 

The double integral of dn s and di?i yields a 8- function of S(p — p' '„), and Eq. (|55|) is simplified 
as: 

V(px,z;p2,z) 

= J dO/(n)G(|pi|,^)G(|p 2 |,^)£ dp A\p )G(\p \,^-) 

dpiG{\M^[-A-\])e-^ + ^ 
c s d J 

W i G(\pM [^[i + i - i]) e-^^K (56) 

lens C S Cl J 

We consider the following two cases: 

Case (I) on the imaging plane and p\= pi= P- 



In this case, Eq. (|56|) is simplified as 



*;?,*) oc / dp o A 2 (p )G(\p \,^) dpte-^^y?' dp^e-^^^'' 

J obj J J 

x{fdQ /(fi) e -« l K&+£MftV,)]} (57) 
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where we have used ui s = uj p /2 + Q and to s = lo p /2 — O following u) s + LUi — uj p . The integral 
of dil gives a 5-function of S[(— + — ){pi — p'i)] while taking the integral to infinity with a 
constant /(f2). This result indicates again that the propagation-diffraction of the signal and 
the idler are not independent. The "two-photon diffraction" couples the two integrals in p Q and 
p' a as well as the two integrals in pi and p'i and thus gives the function 



Po 



GW(p,P)* / dp A 2 (p )e^^ - VS ° C (58) 



obj 



( R u p 

\s C 



Po + ^ 



which indicates that a coherent image (see Eq. (|52[)) magnified by a factor of m = Si/s a is 
reproduced on the image plane by joint-detection or by two-photon absorption. 

In Eq. (|58[) . the point-spread function is characterized by the pump wavelength X p — X si /2: 
hence, the point-spread function is half the width of the (first order) classical case (Eqs. (j52")l and 
(|5ip ). An entangled two-photon state thus gives an image in joint-detection with double spatial 
resolution when compared to the image obtained in classical imaging. Moreover, the spatial 
resolution of the two-photon image obtained by perfect SPDC radiation is further improved 
because it is determined by the function 2J±(x) / 'x 2 , which is much narrower than the somb(x). 

It is interesting to see that, different from the classical case, the frequency integral over 
Aw s ~ oo does not give any blurring, but rather enhances the spatial resolution of the two- 
photon image. 

Case (II): on the Fourier transform plane and p\ = pi= P- 

The detectors are now placed in the focal plane, i.e., d = f. In this case, the spatial effective 
two-photon wavefunction ^(p, z; p, z) becomes: 

*(p>;pU) oc fdSlf(Q) I dp A 2 (p )G(\p \,^) I dp1G(| P 1|A)e-^fe + f^ 

J J obj cs o J lens cs o 

x f dp^Gdp'il^ )e^fe + f>A (59) 

J lens CS 

We will first evaluate the two integrals over the lens. To simplify the mathematics we approxi- 
mate the integral to infinity. Differing from the calculation for imaging resolution, the purpose 
of this evaluation is to determine the Fourier transform. Thus, the approximation of an infinite 
lens is appropriate. By applying Eq. 3D , the two integrals over the lens contribute the 
following function of p to the integral of dp in Eq. (|59l) : 

CG{\p \,-^)e-^°-? 

CS 

where C absorbs all constants including a phase factor G(\p\,— c fi% )■ Replacing the two 
integrals of dpi and dp\ in Eq. ([55]) with this result, we obtain: 

V(p,z;p,z) cx fdflf(Q) [ dp A 2 (p )e~^^°KT niJ] {A 2 (p )}, (60) 

J J obj at 

which is the Fourier transform of the object-aperture function. When the two photodetectors 
scan together (i.e., pi = p2 = p), the second-order transverse correlation G^ 2 ' (p, z; p, z) , where 
z = s + f, is reduced to: 

GW(p,z;p,z) cx \^ pl {A 2 (p )}\ 2 . (61) 

Thus, by replacing classical light with entangled two-photon sources, in the double-slit setup of 
Fig. El a Young's double-slit interference/diffraction pattern with twice the interference modu- 
lation and half the pattern width, compared to that of classical light at wavelength X s i = 2A p , 
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Figure 7: (a) Two-photon Fourier transform of a double-slit. The light source was a collinear 
degenerate SPDC of A S; i = 916nm. (b) Classical Fourier transform of the same double-slit. 
A classical light source of A = 916nm was used. 



is observed in the joint detection. This effect has also been examined in a recent "quantum 
lithography" experiment [20] . 

Due to the lack of two-photon sensitive material, the first experimental demonstration of 
quantum lithography was measured on the Fourier transform plane, instead of the image plane. 
Two point-like photon counting detectors were scanned jointly, similar to the setup illustrated 
in Fig. [H for the observation of the interference/diffraction pattern of Eq. (|61|) . The published 
experimental result is shown in Fig. [7] [20] . It is clear that the two-photon Young's double- 
slit interference-diffraction pattern has half the width with twice the interference modulation 
compared to that of the classical case although the wavelengths are both 916nm. 

Following linear Fourier optics, it is not difficult to see that, with the help of another lens 
(equivalently building a microscope), one can transform the Fourier transform function of the 
double-slit back onto its image plane to observe its image with twice the spatial resolution. 

The key to understanding the physics of this experiment is again through entangled nature 
of the signal-idler two-photon system. As we have discussed earlier, the pair is always emitted 
from the same point on the output plane of the source, thus always passing the same slit 
together if the double-slit is placed close to the surface of the nonlinear crystal. There is no 
chance for the signal-idler pair to pass different slits in this setup. In other words, each point 
of the object is "illuminated" by the pair "together" and the pair "stops" on the image plane 
"together" . The point- "spot" correspondence between the object and image planes are based on 
the physics of two-photon diffraction, resulting in a twice narrower Fourier transform function 
in the Fourier transform plane and twice the image resolution in the image plane. The unfolded 
schematic setup, which is shown in Fig. [8j may be helpful for understanding the physics. It is 
not difficult to calculate the interference-diffraction function under the experimental condition 
indicated in Fig. O The non-classical observation is due to the superposition of the two-photon 
amplitudes, which are indicated by the straight lines connecting D\ and Di- The two-photon 
diffraction, which restricts the spatial resolution of a two-photon image, is very different from 
that of classical light. Thus, there should be no surprise in having an improved spatial resolution 
even beyond the classical limit. 
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Figure 8: Unfolded experimental setup. The joint measurement is on the Fourier transform 
plane. Each point of the object is "illuminated" by the signal-idler pair "together" , result- 
ing in twice narrower interference-diffraction pattern width in the Fourier transform plane 
through the joint detection of the signal-idler pair, equivalent to the use of classical light of 
A/2. 



It is worthwhile to emphasize the following important aspects of physics in this simplified 
illustration: 

(1) The goal of lithography is the reproduction of demagnified images of complicated patterns. 
The sub-wavelength interference feature does not necessarily translate into an improvement of 
the lithographic performance. In fact, the Fourier transform argument works for imaging setups 
only; sub- wavelength interference in a Mach-Zchnder type interferometer, for instance, does not 
necessarily lead to an image. 

(2) In the imaging setup, it is the peculiar nature of the entangled N-photon system that allows 
one to generate an image with N-times the spatial resolution: the entangled photons come out 
from one point of the object plane, undergo N-photon diffraction, and stop in the image plane 
within a N-times narrower spot than that of classical imaging. The historical experiment by 
D'Angelo et al, in which the working principle of quantum lithography was first demonstrated, 
has taken advantage of the entangled two-photon state of SPDC: the signal-idler photon pair 
comes out from either the upper slit or the lower slit that is in the object plane, undergoes 
two-photon diffraction, and stops in the image plane within a twice narrower image than that 
of the classical one. It is easy to show that a second Fourier transform, by means of the use 
of a second lens to set up a simple microscope, will produce an image on the image plane with 
double spatial resolution. 

(3) Certain "clever" tricks allow the production of doubly modulated interference patterns by 
using classical light in joint photo-detection. These tricks, however, may never be helpful for 
imaging. Thus, they may never be useful for lithography. 
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5 Ghost imaging 



The nonlocal position-position and momentum-momentum EPR correlation of the entan- 
gled two-photon state of SPDC was successfully demonstrated in 1995 [21] inspired by the 
theory of Klyshko [22] The experiment was immediately named as "ghost imaging" in the 
physics community due to its surprising nonlocal nature. The important physics demonstrated 
in the experiment, however, may not be the so called "ghost". Indeed, the original purpose 
of the experiment was to study the EPR correlation in position and in momentum and to 
test the EPR inequality of Eq. ^ for the entangled signal- idler photon pair of SPDC [TS] [13] • 
The experiments of "ghost imaging" [2Tj and "ghost interference" [24] together stimulated the 
foundation of quantum imaging in terms of geometrical and physical optics. 
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Figure 9: Schematic set-up of the "ghost" image experiment. 



The schematic setup of the "ghost" imaging experiment is shown in Fig. [9] A CW laser is 
used to pump a nonlinear crystal, which is cut for degenerate type-II phase matching to produce 
a pair of orthogonally polarized signal (e-ray of the crystal) and idler (o-ray of the crystal) 
photons. The pair emerges from the crystal as collinear, with lo s = cjj = uj p /2. The pump is 
then separated from the signal-idler pair by a dispersion prism, and the remaining signal and 
idler beams are sent in different directions by a polarization beam splitting Thompson prism. 
The signal beam passes through a convex lens with a 400mm focal length and illuminates a 
chosen aperture (mask). As an example, one of the demonstrations used the letters "UMBC" 
for the object mask. Behind the aperture is the "bucket" detector package D\, which consists 
of a short focal length collection lens in whose focal spot is an avalanche photodiode. D\ is 
mounted in a fixed position during the experiment. The idler beam is met by detector package 
D2, which consists of an optical fiber whose output is mated with another avalanche photodiode. 
The input tip of the fiber is scanned in the transverse plane by two step motors. The output 
pulses of each detector, which are operating in photon counting mode, are sent to a coincidence 
counting circuit for the signal-idler joint detection. 

By recording the coincidence counts as a function of the fiber tip's transverse plane coordi- 
nates, the image of the chosen aperture (for example, "UMBC") is observed, as reported in Fig. 
[TU1 It is interesting to note that while the size of the "UMBC" aperture inserted in the signal 
beam is only about 3.5mm x 7mm, the observed image measures 7mm x 14mm. The image 
is therefore magnified by a factor of 2. The observation also confirms that the focal length of 
the imaging lens, /, the aperture's optical distance from the lens, S Q1 and the image's optical 
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Figure 10: (a) A reproduction of the actual aperature "UMBC" placed in the signal beam, 
(b) The image of "UMBC" : coincidence counts as a function of the fiber tip's transverse 
plane coordinates. The step size is 0.25mm. The image shown is a "slice" at the half 
maximum value. 

distance from the lens, Si (which is from the imaging lens going backward along the signal 
photon path to the two-photon source of the SPDC crystal then going forward along the path 
of idler photon to the image), satisfy the Gaussian thin lens equation. In this experiment, S 
was chosen to be S = 600mm, and the twice magnified clear image was found when the fiber 
tip was on the plane of Si — 1200mm. While D 2 was scanned on other transverse planes not 
defined by the Gaussian thin lens equation, the images blurred out. 

The measurement of the signal and the idler subsystem themselves are very different. The 
single photon counting rate of D 2 was recorded during the scanning of the image and was found 
fairly constant in the entire region of the image. This means that the transverse coordinate 
uncertainty of either signal or idler is considerably large compared to that of the transverse 
correlation of the entangled signal-idler photon pair: Axi (Ayi) and Ax 2 (Aj/2) are much 
greater than A(xi - x 2 ) (A(yi - y 2 )). 

The EPR (^-functions, S(p s — pi) and S(k s + k{) in transverse dimension, are the key to 
understanding this interesting phenomenon. In degenerate SPDC, although the signal- idler 
photon pair has equal probability to be emitted from any point on the output surface of the 
nonlinear crystal, the transverse position (5-function indicates that if one of them is observed 
at one position, the other one must be found at the same position. In other words, the pair 
is always emitted from the same point on the output plane of the two-photon source. The 
transverse momentum (5-function, defines the angular correlation of the signal-idler pair: the 
transverse momenta of a signal-idler amplitude are equal but pointed in opposite directions: 
In other words, the two-photon amplitudes are always existing at roughly equal 
yet opposite angles relative to the pump. This then allows for a simple explanation of the 
experiment in terms of "usual" geometrical optics in the following manner: we envision the 
nonlinear crystal as a "hinge point" and "unfold" the schematic of Fig. [9] into that shown in 
Fig. 1111 The signal-idler two-photon amplitudes can then be represented by straight lines (but 
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Figure 11: An unfolded setup of the "ghost" imaging experiment, which is helpful for un- 
derstanding the physics. Since the two-photon "light" propagates along "straight-lines" , it 
is not difficult to find that any geometrical light point on the subject plane corresponds to 
an unique geometrical light point on the image plane. Thus, a "ghost" image of the subject 
is made nonlocally in the image plane. Although the placement of the lens, the object, and 
detector D2 obeys the Gaussian thin lens equation, it is important to remember that the 
geometric rays in the figure actually represent the two-photon amplitudes of an entangled 
photon pair. The point to point correspondence is the result of the superposition of these 
two-photon amplitudes. 



keep in mind the different propagation directions) and therefore, the image is well produced 
in coincidences when the aperture, lens, and fiber tip are located according to the Gaussian 
thin lens equation of Eq. ([5]) . The image is exactly the same as one would observe on a screen 
placed at the fiber tip if detector D\ were replaced by a point-like light source and the nonlinear 
crystal by a reflecting mirror. 

Following a similar analysis in geometric optics, it is not difficult to find that any geometri- 
cal "light spot" on the subject plane, which is the intersection point of all possible two-photon 
amplitudes coming from the two-photon light source, corresponds to a unique geometrical "light 
spot" on the image plane, which is another intersection point of all the possible two-photon 
amplitudes. This point to point correspondence made the "ghost" image of the subject-aperture 
possible. Despite the completely different physics from classical geometrical optics, the remark- 
able feature is that the relationship between the focal length of the lens, /, the aperture's 
optical distance from the lens, S , and the image's optical distance from the lens, Si, satisfy 
the Gaussian thin lens equation: 

1 1 _ 1 

So Si f 

Although the placement of the lens, the object, and the detector D 2 obeys the Gaussian thin lens 
equation, it is important to remember that the geometric rays in the figure actually represent 
the two-photon amplitudes of a signal-idler photon pair and the point to point correspondence 
is the result of the superposition of these two-photon amplitudes. The "ghost" image is a 
realization of the 1935 EPR gedankenexperiment. 

Now we calculate G^ 2 '{p 01 pi) for the "ghost" imaging experiment, where p and pi are the 
transverse coordinates on the object plane and the image plane. We will show that there exists 
a ^-function like point-to-point relationship between the object plane and the image plane, i.e., 
if one measures the signal photon at a position of p on the object plane the idler photon can be 
found only at a certain unique position of pi on the image plane satisfying 5(mp — pi), where 
m = — (si/s ) is the image-object magnification factor. After demonstrating the (5-function, we 
show how the object-aperture function of A(p Q ) is transfered to the image plane as a magnified 
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image A(pi/m). Before showing the calculation, it is worthwhile to emphasize again that the 
"straight lines" in Fig. [TT] schematically represent the two-photon amplitudes belonging to a 
pair of signal-idler photon. A "click-click" joint measurement at (ri, ti), which is on the object 
plane, and (T2, £2), which is on the image plane, in the form of an EPR <5-function, is the result 
of the coherent superposition of all these two-photon amplitudes. 




Figure 12: In arm-1, the signal propagates freely over a distance d\ from the output plane 
of the source to the imaging lens, then passes an object aperture at distance s a , and then is 
focused onto photon counting detector D\ by a collection lens. In arm-2, the idler propagates 
freely over a distance cfo from the output plane of the source to a point-like photon counting 
detector Z?2- 

We follow the unfolded experimental setup shown in Fig. [T2] to establish the Green's func- 
tions q(k s , lo Si p , z ) and g(Kj, Wj, P2, 22). In arm-1, the signal propagates freely over a distance 
di from the output plane of the source to the imaging lens, then passes an object aperture at 
distance s , and then is focused onto photon counting detector D\ by a collection lens. We 
will evaluate g(K s ,Lu s , p Q , z ) by propagating the field from the output plane of the two-photon 
source to the object plane. In arm-2, the idler propagates freely over a distance c?2 from the 
output plane of the two-photon source to a point-like detector D2. g(Ki,uii, P2, 22) is thus a free 
propagator. 

(I) Arm-1 (source to object): 

The optical transfer function or Green's function in arm-1, which propagates the field from 
the source plane to the object plane, is given by: 

g{K s ,u s ;p ,z = <fi + s a ) 

= e^*° [ dp t ( dpsl^e^GQps-pil^-)} 

Jlens J A L 2-Kcdx cd X > 

x{G(|,1|,^)}{-^lG(|p1-p |,^)}, (62) 
I CJ J I 2ttcs cs J 

where ps and pi are the transverse vectors defined, respectively, on the output plane of the 
source and on the plane of the imaging lens. The terms in the first and third curly brackets 
in Eq. (|62[) describe free space propagation from the output plane of the source to the imaging 
lens and from the imaging lens to the object plane, respectively. The function G(\pi\, ^j) in the 
second curly brackets is the transformation function of the imaging lens. Here, we treat it as a 
thin-lens: G(\pi\, = e^W l^l 2 . 
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(II) Arm-2 (from source to image): 

In arm-2, the idler propagates freely from the source to the plane of D 2 , which is also the 
plane of the image. The Green's function is thus: 

g (H i ,iv i ;p 2 ,z 2 = d 2 ) = ^-e i ^ f d p' s G{\p' s -p 2 \^)e^s (63) 
2ncd 2 J a c «2 

where p' s and p 2 are the transverse vectors defined, respectively, on the output plane of the 
source, and on the plane of the photo-dector D 2 . 

(III) fy(p ,pi) (object plane - image plane): 

To simplify the calculation and to focus on the transverse correlation, in the following 
calculation we assume degenerate (uj s = uJi = us) and collinear SPDC. The transverse two- 
photon effective wavefunction ^{p , P2) is then evaluated by substituting the Green's functions 
g(K s ,uj; p , z ) and g(Ki, to; p 2 , z 2 ) into the expression given in Eq. (f3"B"|): 

oc / dK s dKiS(K s + Ki) g(K s ,Lu; p , z ) g(Ki,oj; p 2l z 2 ) 
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x G(\pi\, -) G{\ px -Pol—) dp> s e^s G(\ p> s -p 2 \, (64) 
cj cs j A ca 2 

where we have ignored all the proportional constants. Completing the double integral of di? s 
and dn s 

J dn s dRi 5{k s + Hi) e l ^P s „ 5(p s - p r s ), (65) 

Eq. l[Mp becomes: 
*(Po,p 2 ) 

oc / dpi dp s G(\p 2 ~ ps\,^)G(\p s - pi\,-^-)G(\pi\,^-)G(\pi- p 
Jlens Ja cd 2 cdi c/ 
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We then apply the properties of the Gaussian functions of Eq. flyl— 3| ) and complete the integral 
on dps by assuming the transverse size of the source is large enough to be treated as infinity. 

*(p ,p 2 )oc/ dpi G(\ p 2 - pi I, — ) G(\pi\, ^j) G(\ pi - p I, — ). (66) 

Although the signal and idler propagate to different directions along two optical arms, Inter- 
estingly, the Green function in Eq. (|66|) is equivalent to that of a classical imaging setup, if 
we imagine the fields start propagating from a point p on the object plane to the lens and 
then stop at point p 2 on the imaging plane. The mathematics is consistent with our previous 
qualitative analysis of the experiment. 

The integral on dpi yields a point-to-point relationship between the object plane and the 
image plane that is defined by the Gaussian thin-lens equation: 

dpi G(\ pi\,-[- + -- h) e -^fe+^ oc S(p + h) (67) 

lens C So Si J TJl 

where the integral is approximated to infinity and the Gaussian thin-lens equation of l/s + 
1/ Si = 1// is applied. We have also defined m = Si/s as the magnification factor of the imaging 
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system. The function S(p + pi/m) indicates that a point p D on the object plane corresponds 
to a unique point pi on the image plane. The two vectors point in opposite directions and the 
magnitudes of the two vectors hold a ratio of in = \pi\/\p \. 

If the finite size of the imaging lens has to be taken into account (finite diameter D), the 
integral yields a point-spread function of somb(x) : 

/ dp ie - l ^ + ^ y ? 1 oc somb( - -\p + &]) (68) 

J lens Vs C m J 

where somb(x) = 2Ji(x)/x, J\{x) is the first-order Bessel function and R/s is named as the 
numerical aperture. The point-spread function turns the point-to-point correspondence between 
the object plane and the image plane into a point-to- "spot" relationship and thus limits the 
spatial resolution. This point has been discussed in detail in the last section. 

Therefore, by imposing the condition of the Gaussian thin-lens equation, the transverse 
two-photon effective wavefunction is approximated as a 8 function 

<xS(p + ^) (69) 
m 

where p and pi, again, are the transverse coordinates on the object plane and the image 
plane, respectively, defined by the Gaussian thin-lens equation. Thus, the second-order spatial 
correlation function G^ (p , p%) turns out to be: 

G^(Po, Pi) = | We Pi) I 2 * I 5{Po + | 2 . (70) 

m 

Eq. ([70)) indicates a point to point EPR correlation between the object plane and the image 
plane, i.e., if one observes the signal photon at a position p on the object plane, the idler photon 
can only be found at a certain unique position pi on the image plane satisfying S(p + Pi/m) 
with m = Si/s a . 



We now include an object-aperture function, a collection lens and a photon counting de- 
tector Di into the optical transfer function of arm-1 as shown in Fig.[5J 

We will first treat the collection-lens-Di package as a "bucket" detector. The "bucket" de- 
tector integrates all ^(p , Pi) which passes the object aperture A{p ) as a joint photo-detection 
event. This process is equivalent to the following convolution : 

Hi, 2 cx / dp |A(p )| 2 |*(p ,p l )| 2 -|^( — )| 2 (71) 



obj 



in 



where, again, D2 is scanning in the image plane, pi= Pi- Eq. J71|) indicates a magnified (or de- 
magnified) image of the object-aperture function by means of the joint-detection events between 
distant photodetectors D\ and D2. The "-" sign in A(—pi/m) indicates opposite orientation of 
the image. The model of the "bucket" detector is a good and realistic approximation. 

Now we consider a detailed evaluation by including the object-aperture function, the col- 
lection lens and the photon counting detector D\ into arm-1. The Green's function of Eq. (|62[) 
becomes: 

g(K s , uj s \Pi, zi =dt + s + f coU ) 
- [ dp [ dpt ( dp s { Z!^. e «i*G(| p s - p t |, 

Job] hens J A L ^Cd X cd X i 

x G{\pi\, ^) [ ^p^-G(\pi - p I, —)XA(p ) 
cj I Zircso cs > 

x G(|p c |, J£-){ — ^- G(\p J£-)\ (72) 

CJcoll ^ tncjcoll CJcoll J 
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where f co u is the focal-length of the collection lens and D\ is placed on the focal point of 
the collection lens. Repeating the previous calculation, we obtain the transverse two-photon 
effective wavefunction: 

*(pi,p 2 )oc/ dp A(p ) S{p + ^) = A(p ) ® 8(p + (73) 
Jobj m m 

where <8> means convolution. Notice, in Eq. (|T3() we have ignored the phase factors which have 
no contribution to the formation of the image. The joint detection counting rate, i?i,2, between 
photon counting detectors D\ and D2 is thus: 

Rx,2 oc G^(p!,p 2 ) oc I A(p a ) ® S(p + ^) I 2 - I Af^) | 2 (74) 
where, again, P2 — (H- 

As we have discussed earlier, the point-to-point EPR correlation is the result of the coherent 
superposition of a special selected set of two-photon states. In principle, one signal-idler pair 
contains all the necessary two-photon amplitudes that generate the ghost image - a nonclassical 
characteristic which we name as a two-photon coherent image. 



6 Popper's experiment 

In quantum mechanics, one can never expect to measure both the precise position and mo- 
mentum of a particle simultaneously. It is prohibited. We say that the quantum observable 
"position" and "momentum" are "complementary" because the precise knowledge of the po- 
sition (momentum) implies that all possible outcomes of measuring the momentum (position) 
are equally probable. 

Karl Popper, being a "metaphysical realist" , however, took a different point of view. In his 
opinion, the quantum formalism could and should be interpreted realistically: a particle must 
have a precise position and momentum [25) . This view was shared by Einstein. In this regard, he 
invented a thought experiment in the early 1930's aimed to support his realistic interpretation 
of quantum mechanics [26j . What Popper intended to show in his thought experiment is that a 
particle can have both precise position and momentum simultaneously through the correlation 
measurement of an entangled two-particle system. 

Similar to EPR's gedankenexperiment, Popper's thought experiment is also based on the 
feature of two-particle entanglement: if the position or momentum of particle 1 is known, the 
corresponding observable of its twin, particle 2, is then 100% determined. Popper's original 
thought experiment is schematically shown in Fig. 1131 A point source S, positronium as Popper 
suggested, is placed at the center of the experimental arrangement from which entangled pairs 
of particles 1 and 2 are emitted in opposite directions along the respective positive and negative 
x-axes towards two screens A and B. There are slits on both screens parallel to the y-axis and 
the slits may be adjusted by varying their widths Ay. Beyond the slits on each side stand an 
array of Geiger counters for the joint measurement of the particle pair as shown in the figure. 
The entangled pair could be emitted to any direction in An solid angles from the point source. 
However, if particle 1 is detected in a certain direction, particle 2 is then known to be in the 
opposite direction due to the momentum conservation of the pair. 

First, let us imagine the case in which slits A and B are both adjusted very narrowly. In this 
circumstance, particle 1 and particle 2 experience diffraction at slit A and slit B, respectively, 
and exhibit greater Ap y for smaller Ay of the slits. There seems to be no disagreement in this 
situation between Copenhagen and Popper. 

Next, suppose we keep slit A very narrow and leave slit B wide open. The main purpose 
of the narrow slit A is to provide the precise knowledge of the position y of particle 1 and 
this subsequently determines the precise position of its twin (particle 2) on side B through 
quantum entanglement. Now, Popper asks, in the absence of the physical interaction with an 
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(a) 




Figure 13: Popper's thought experiment. An entangled pair of particles are emitted from 
a point source with momentum conservation. A narrow slit on screen A is placed in the 
path of particle 1 to provide the precise knowledge of its position on the j/-axis and this also 
determines the precise y-position of its twin, particle 2, on screen B. (a) Slits A and B are 
both adjusted very narrowly, (b) Slit A is kept very narrow and slit B is left wide open. 

actual slit, does particle 2 experience a greater uncertainty in Ap y due to the precise knowledge 
of its position? Based on his beliefs, Popper provides a straightforward prediction: particle 2 
must not experience a greater Ap y unless a real physical narrow slit B is applied. However, 
if Popper's conjecture is correct, this would imply the product of Ay and Ap y of particle 2 
could be smaller than h (Ay Ap y < h). This may pose a serious difficulty for Copenhagen and 
perhaps for many of us. On the other hand, if particle 2 going to the right does scatter like its 
twin, which has passed though slit A, while slit B is wide open, we are then confronted with 
an apparent action- at- a- distance! 

The use of a "point source" in Popper's proposal has been criticized historically as the 
fundamental mistake Popper made [27) . It is true that a point source can never produce a pair 
of entangled particles which preserves the EPR correlation in momentum as Popper expected. 
However, notice that a "point source" is not a necessary requirement for Popper's experiment. 
What is required is a precise position-position EPR correlation: if the position of particle 1 is 
precisely known, the position of particle 2 is 100% determined. As we have shown in the last 
section, "ghost" imaging is a perfect tool to achieve this. 

In 1998, Popper's experiment was realized with the help of two-photon "ghost" imaging 
[28] . Fig. [14] is a schematic diagram that is useful for comparison with the original Popper's 
thought experiment. It is easy to see that this is a typical "ghost" imaging experimental setup. 
An entangled photon pair is used to image slit A onto the distant image plane of "screen" B. In 
the setup, s is chosen to be twice the focal length of the imaging lens LS, s — 2/. According 
to the Gaussian thin lens equation, an equal size "ghost" image of slit A appears on the two- 
photon image plane at Sj = 2/. The use of slit A provides a precise knowledge of the position 
of photon 1 on the j/-axis and also determines the precise y-position of its twin, photon 2, on 
screen B by means of the two-photon "ghost" imaging. The experimental condition specified in 
Popper's experiment is then achieved. When slit A is adjusted to a certain narrow width and 
slit B is wide open, slit A provides precise knowledge about the position of photon 1 on the 
y-axis up to an accuracy Ay which equals the width of slit A, and the corresponding "ghost 
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Figure 14: Modified version of Popper's experiment. An entangled photon pair is generated 
by SPDC. A lens and a narrow slit A are placed in the path of photon 1 to provide the 
precise knowledge of its position on the y-axis and also to determine the precise y-position 
of its twin, photon 2, on screen B by means of two-photon "ghost" imaging. Photon counting 
detectors D\ and D2 are used to scan in y-directions for joint detections, (a) Slits A and B 
are both adjusted very narrowly, (b) Slit A is kept very narrow and slit B is left wide open. 



image" of pinhole A at screen B determines the precise position y of photon 2 to within the 
same accuracy Ay. Ap y of photon 2 can be independently studied by measuring the width of 
its "diffraction pattern" at a certain distance from "screen" B. This is obtained by recording 
coincidences between detectors D\ and D2 while scanning detector D2 along its y-axis, which 
is behind screen B at a certain distance. 

Figure [TBI is a conceptual diagram to connect the modified Popper's experiment with two- 
photon "ghost" imaging. In this unfolded "ghost" imaging setup, we assume the entangled 
signal- idler photon pair holds a perfect transverse momentum correlation with k s + ki ~ 0, 
which can be easily realized in SPDC. In this experiment, we have chosen s D = Sj = 2/. Thus, 
an equal size "ghost" image of slit A is expected to appear on the image plane of screen B. 

The detailed experimental setup is shown in Fig |16l with indications of the various distances. 
A CW Argon ion laser line of X p — 351. Inm is used to pump a 3mm long beta barium borate 
(BBO) crystal for type-II SPDC to generate an orthogonally polarized signal-idler photon pair. 
The laser beam is about 3mm in diameter with a diffraction limited divergence. It is important 
to keep the pump beam a large size so that the transverse phase-matching condition, k s + ki ~ 
(k p = 0), is well reinforced in the SPDC process, where kj (j = s, i) is the transverse wavevector 
of the signal (s) and idler (i), respectively. The collinear signal-idler beams, with A s = = 
702. 2nm — 2X p are separated from the pump beam by a fused quartz dispersion prism, and 
then split by a polarization beam splitter PBS. The signal beam (photon 1) passes through the 
converging lens LS with a 500mm focal length and a 25mm diameter. A 0.16mm slit is placed 
at location A which is 1000mm (= 2/) behind the lens LS. A short focal length lens is used 
with Di for focusing the signal beam that passes through slit A. The point-like photon counting 
detector D2 is located 500mm behind "screen B" . "Screen B" is the image plane defined by the 
Gaussian thin lens equation. Slit B, either adjusted as the same size as that of slit A or opened 
completely, is placed to coincide with the "ghost" image. The output pulses from the detectors 
are sent to a coincidence circuit. During the measurements, detector D\ is fixed behind slit A 
while detector D2 is scanned on the y-axis by a step motor. 
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Figure 15: An unfolded schematic of ghost imaging. We assume the entangled signal-idler 
photon pair holds a perfect momentum correlation #(k s + kj) ~ 0. The locations of the slit 
A, the imaging lens LS, and the "ghost" image must be governed by the Gaussian thin lens 
equation. In this experiment, we have chosen s = Sj = 2/. Thus, the "ghost" image of slit 
A is expected to be the same size as that of slit A. 

Measurement 1: Measurement 1 studied the case in which both slits A and B were 
adjusted to be 0.16mm. The y-coordinate of D\ was chosen to be (center) while D 2 was 
allowed to scan along its y-axis. The circled dot data points in Fig. [17] show the coincidence 
counting rates against the y-coordinates of D 2 . It is a typical single-slit diffraction pattern with 
Ay Ap y — h. Nothing is special in this measurement except that we have learned the width 
of the diffraction pattern for the 0.16mm slit and this represents the minimum uncertainty 
of Ap y . We should emphasize at this point that the single detector counting rate of D2 as a 
function of its position y is basically the same as that of the coincidence counts except for a 
higher counting rate. 

Measurement 2: The same experimental conditions were maintained except that slit B 
was left wide open. This measurement is a test of Popper's prediction. The y-coordinate of D\ 
was chosen to be (center) while D 2 was allowed to scan along its y-axis. Due to the entangled 
nature of the signal- idler photon pair and the use of a coincidence measurement circuit, only 
those twins which have passed through slit A and the "ghost image" of slit A at screen B with 
an uncertainty of Ay = 0.16mm (which is the same width as the real slit B we have used in 
measurement 1) would contribute to the coincidence counts through the joint detection of D\ 
and D 2 . The diamond dot data points in Fig. [T7| report the measured coincidence counting 
rates against the y coordinates of D 2 . The measured width of the pattern is narrower than 
that of the diffraction pattern shown in measurement 1. It is also interesting to notice that the 
single detector counting rate of D 2 keeps constant in the entire scanning range, which is very 
different from that in measurement 1. The experimental data has provided a clear indication 
of Ay Ap y < h in the joint measurements of the entangled photon pairs. 

Given that Ay Ap y < h, is this a violation of the uncertainty principle? Does quantum 
mechanics agree with this peculiar experimental result? If quantum mechanics does provide a 
solution with Ay Ap y < h for photon 2. We would indeed be forced to face a paradox as EPR 
had pointed out in 1935. 

Quantum mechanics does provide a solution that agrees with the experimental result. How- 
ever, the solution is for a joint measurement of an entangled photon pair that involves both 
photon 1 and photon 2, but not just for photon 2 itself . 

We now examine the experimental results with the quantum mechanical calculation by 
adopting the formalisms from the ghost image experiment with two modifications: 

Case (I): slits A = 0.16mm, slit B = 0.16mm. 
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Figure 16: Schematic of the experimental setup. The laser beam is about 3mm in diameter. 
The "phase- matching condition" is well reinforced. Slit A (0.16mm) is placed 1000mm = 2/ 
behind the converging lens, LS (/ = 500mm). The one-to-one "ghost image" (0.16mm) of 
slit A is located at B. The optical distance from LS in the signal beam taken as back 
through PBS to the SPDC crystal (£>i = 255mm) and then along the idler beam to "screen 
B" (&2 = 745mm) is 1000mm = 2/ (b = 6i + b 2 ). 

This is the experimental condition for measurement one: slit B is adjusted to be the same as 
slit A. There is nothing surprising about this measurement. The measurement simply provides 
us with the knowledge for Ap y of photon 2 caused by the diffraction of slit B (Ay = 0.16mm). 
The experimental data shown in Fig. [T7| agrees with the calculation. Notice that slit B is about 
745mm away from the 3mm two-photon source, the angular size of the light source is roughly 
the same as A/Ay, AO ~ A/Ay, where A = 702nm is the wavelength and Ay = 0.16mm is 
the width of the slit. The calculated diffraction pattern is very close to that of the "far-field" 
Fraunhofer diffraction of a 0.16mm single-slit. 

Case (II): slit A = 0.16mm, slits B ~ oo (wide open). 

Now we remove slit B from the ghost image plane. The calculation of the transverse effective 
two-photon wavefunction and the second-order correlation is the same as that of the ghost image 
except the observation plane of D 2 is moved behind the image plane to a distance of 500mm. 
The two-photon image of slit A is located at a distance s, = 2/ = 1000mm (b\ + 62) from the 
imaging lens, in this measurement D 2 is placed at d = 1500mm from the imaging lens. The 
measured pattern is simply a "blurred" two-photon image of slit A. The "blurred" two-photon 
image can be calculated from Eq. ((75)) which is a slightly modified version of Eq. (|66[) 



where d is the distance between the imaging lens and D 2 . In this measurement, D 2 was placed 
500mm behind the image plane, i.e., d = s$ + 500mm. The numerically calculated "blurred" 
image, which is narrower then that of the diffraction pattern of the 0.16mm slit B, agrees with 
the measured result of Fig. [T7] within experimental error. 

The measurement does show a result of Ay Ap y < h. The measurement, however, has 
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Figure 17: The observed coincidence patterns. The y-coordinate of D\ was chosen to be 
(center) while Di was allowed to scan along its y-axis. Circled dot points: Slit A = Slit B = 
0.16mm. Diamond dot points: Slit A = 0.16mm, Slit B wide open. The width of the sine 
function curve fitted by the circled dot points is a measure of the minimum Ap y diffracted 
by a 0.16mm slit. 



nothing to do with the uncertainty relation, which governs the behavior of photon 2 (the idler) . 
Popper and EPR were correct in the prediction of the outcomes of their experiments. Popper 
and EPR, on the other hand, made the same error by applying the results of two-particle physics 
to the explanation of the behavior of an individual subsystem. 

In both the Popper and EPR experiments, the measurements are "joint detection" between 
two detectors applied to entangled states. Quantum mechanically, an entangled two-particle 
state only provides the precise knowledge of the correlations of the pair. The behavior of 
"photon 2" observed in the joint measurement is conditioned upon the measurement of its 
twin. A quantum must obey the uncertainty principle but the "conditional behavior" of a 
quantum in an entangled two-particle system is different in principle. We believe paradoxes are 
unavoidable if one insists the conditional behavior of a particle is the behavior of the particle. 
This is the central problem in the rationale behind both Popper and EPR. Ay Ap y > h is not 
applicable to the conditional behavior of either "photon 1" or "photon 2" in the cases of Popper 
and EPR. 

The behavior of photon 2 being conditioned upon the measurement of photon 1 is well 
represented by the two-photon amplitudes. Each of the straight lines in the above discus- 
sion corresponds to a two-photon amplitude. Quantum mechanically, the superposition of 
these two-photon amplitudes are responsible for a "click-click" measurement of the entangled 
pair. A "click-click" joint measurement of the two-particle entangled state projects out certain 
two-particle amplitudes, and only these two-particle amplitudes are featured in the quantum 
formalism. In the above analysis we never consider "photon 1" or "photon 2" individually. 
Popper's question about the momentum uncertainty of photon 2 is then inappropriate. 

Once again, the demonstration of Popper's experiment calls our attention to the important 
message: the physics of an entangled two-particle system must be inherently very different from 
that of individual particles. 



7 Subsystem in an entangled two-photon state 

The entangled EPR two-particle state is a pure state with zero entropy. The precise corre- 
lation of the subsystems is completely described by the state. The measurement, however, is 



36 



not necessarily always on the two-photon system. It is an experimental choice to study a single 
subsystem and to ignore the other. What can be learn about a subsystem from these kinds 
of measurements? Mathematically, it is easy to show that by taking a partial trace of a two- 
particle pure state, the state of each subsystem is in a mixed state with entropy greater than 
zero. One can only learn statistical properties of the subsystems in this kind of measurement. 

In the following, again, we use the signal-idler pair of SPDC as an example to study the 
physics of a subsystem. The two-photon state of SPDC is a pure state that satisfies 

p 2 = p, p=\*)(*\ 

where p is the density operator corresponding to the two-photon state of SPDC. The single 
photon states of the signal and idler 

p s = tri |*) (*| , ft = tr.|*)(¥| 

are not pure states. To calculate the signal (idler) state from the two-photon state, we take a 
partial trace, as usual, summing over the idler (signal) modes. 

We assume a type II SPDC. The orthogonally polarized signal and idler are degenerate 
in frequency around = w s ° = u> p /2. To simplify the discussion, by assuming appropriate 
experimental conditions, we trivialize the transverse part of the state and write the two-photon 
state in the following simplified form: 

|*) = * J dQ $(DLfi) al(u° + fi) ot(w? - fl) |0) 

where $(DLJ7) is a sinc-like function: 

i _ p — tDLfi 
$(DLO) = — 

which is a function of the crystal length L, and the difference of inverse group velocities of the 
signal (ordinary) and the idler (extraordinary), D = l/u — l/u e . The constant is calculated 
from the normalization tr p = | ^) = 1. It is easy to calculate and to find p 2 = p for the 
two-photon state of the signal-idler pair. 

Summing over the idler modes, the density matrix of signal is given by 

p s = ^l J dn |$(fi)| 2 aj(w° + n) |0) (0| a s (cj° s + Q) (76) 

with 

|^)| 2 =sinc 2 ^ 

where all constants coming from the integral have been absorbed into \Po - First, we find 
immediately that p 2 ^ p s . It means the state of the signal is a mixed state (as is the idler). 
Second, it is very interesting to find that the spectrum of the signal depends on the group 
velocity of the idler. This, however, should not come as a surprise, because the state of the 
signal photon is calculated from the two-photon state by summing over the idler modes. 

The spectrum of the signal and idler has been experimentally verified by Strekalov et al 
using a Michelson interferometer in a standard Fourier spectroscopy type measurement [29] . The 
measured interference pattern is shown in Fig. 1181 The envelope of the sinusoidal modulations 
(in segments) is fitted very well by two "notch" functions (upper and lower part of the envelope). 
The experimental data agrees with the theoretical analysis of the experiment. 

The following is a simple calculation to explain the observed "notch" function. We first 
define the field operators: 

EW(t, z d ) = EW(t -^,zo) + (t -^,z ) 
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Figure 18: Experimental data indicated a "double notch" envelope. Each of the doted single 
vertical line contines many cycles of sinusoidal modulation. 



where z d is the position of the photo-detector, zq is the input point of the interferometer, 
t\ = t — — and < 2 = t — ^r, respectively, are the early times before propagating to the photo- 
detector at time t with time delays of Z\jc and 22/c, where Z\ and z 2 are the optical paths in 
arm 1 and arm 2 of the interferometer. We have defined a very general field operator which 
is a superposition of two early fields propagated individually through arm 1 and arm 2 of any 
type of interferometer. The counting rate of the photon counting detector is thus 



tr [p s E^(t,z d )E^(t,z d )] 



dn |$(fi)| 2 \ (0\E (+ \t, z d ) al{uj s +0) |0) 



dn mn)\ 2 1 <o| [£?<+>(* - ^, Zo ) + - z -f, zo)] aM + n) |o) 



oc 1+Re 



dVt sine 



--iOr 



(77) 



where r = (21 — Z2)/c. The Fourier transform of sinc 2 (DLf2/2) has a "notch" shape. It is 
noticed that the base of the "notch" function is determined by parameter DL of the SPDC, 
which is easily confirmed from the experiment. 



Now we turn to another interesting aspect of physics, namely the physics of entropy. In 
classical information theory, the concept of entropy, named as Von Neuman entropy, is defined 
by [30] 

S = -tr(p log p) (78) 

where p is the density operator. It is easy to find that the entropy of the entangled two-photon 
pure state is zero. The entropy of its subsystems, however, are both greater than zero. The value 
of the Von Neuman entropy can be numerically evaluated from the measured spectrum. Note 
that the density operator of the subsystem is diagonal. Taking its trace is simply performing an 
integral over the frequency spectrum with the measured spectrum function. It is straightforward 
to find the entropy of the subsystems S s > 0. This is an expected result due to the statistical 
mixture nature of the subsystem. Considering that the entropy of the two-photon system is zero 
and the entropy of the subsystems are both greater than zero, does this mean that negative 
entropy is present somewhere in the entangled two-photon system? According to classical 
"information theory", for the entangled two-photon system, S s + S 3 \i = 0, where S s ^ is the 
conditional entropy. It is this conditional entropy that must be negative, which means that 
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given the result of a measurement over one particle, the result of a measurement over the other 
must yield negative information [31 j . This paradoxical statement is similar and, in fact, closely 
related to the EPR "paradox" . It comes from the same philosophy as that of the EPR. 

Summary 

The physics of an entangled system is very different from that of either classically inde- 
pendent or correlated systems. We use 2 ^ 1 + 1 to emphasize the nonclassical behavior of an 
entangled two-particle system. The entangled system is characterized by the properties of an 
entangled state which does not specify the state of an individual system, but rather describes 
the correlation between the subsystems. An entangled two particle state is a pure state which 
involves the superposition of a set of "selected" two-particle states, or two-particle quantum 
mechanical amplitudes. Here, the term "selection" stems from the physical laws which govern 
the creation of the subsystems in the source, such as energy or momentum conservation. In- 
terestingly, quantum mechanics allows for the superposition of these local two-particle states 
which have been observed in nature. However, the most surprising physics arises from the joint 
measurement of the two particles when they are released form the source and propagated a large 
distance apart. The two well separated interaction- free particles do not lose their entangled 
properties, i.e., they maintain their "selected" set of two-particle superposition. In this sense 
quantum mechanics allows for the two-particle superposition of well separated particles which 
has, remarkably, also been observed to exist in nature. 

The two-photon state of SPDC is a good example. The nonlinear interaction of sponta- 
neous parametric down-conversion coherently creates a set of mode in pairs that satisfy the 
phase matching conditions of Eq. which is also characteristic of energy and momentum 
conservation. The signal-idler photon pair can be excited to any or all of these coupled modes 
simultaneously, resulting in a superposition of these coupled modes inside of the nonlinear crys- 
tal. The physics behind the two-photon superposition becomes even more interesting when 
the signal-idler pair is separated and propagated a large distance apart outside the nonlinear 
crystal, either through free propagation or guided by optical components. Remarkably the en- 
tangled pair does not lose its entangled properties once the subsystems are interaction free. As 
a result the properties of the entangled two-photon system, such as the EPR correlation or the 
EPR inequalities, are still observable in the joint detection counting rate of the pair, regardless 
of the distance between the two photons as well as the two individual photo-detection events. 
In this situation the superposition of the two-photon amplitudes, corresponding to different yet 
indistinguishable alternative ways of triggering a joint photo-electron event at any distance can 
be regarded as nonlocal. There is no counterpart to such a concept in classical theory and this 
behavior may never be understood in any classical sense. It is with this intent that we use 
2 =/= 1 + 1 to emphasize that the physics of a two-photon is not the same as that of two photons. 



A statement from the author 

This article was originally prepared as lecture notes for my students a few years ago. It was 
also used in 2006 for a conference. My colleagues, friends and students have urged me to include 
it in this archive. They believe that this article is helpful for the general physics and engineering 
community. Truthfully, I have been hesitant because I cannot forget my terrible experience in 
1996 as a result of Pittman's experiment: "Can tow-photon interference be considered the 
interference of two photons?" [21] My email account was bombarded for months. Of course, I 
was happy to have scientific discussions on the subject, but certain types of messages caused 
headaches. For example, an individual attempted to force my laboratory to pay a visit for a 
face-to- face condemnation on my guilt for saying 1 + 1^2. (I truly believe what I said was 
2^1 + 1 and anyone would be able to see the difference by reading this article). Another 
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individual expressed their interests in a law suit because we did not acknowledge that they 
were the first to show "Dirac was mistaken" . (I am definitely sure that we have nothing to 
do with their "discovery". What we said was "Dirac was correct".) I decided to keep quiet. I 
understood that it takes time for people to recognize the truth. 

I have to break my silence now, because we are experiencing the same problem again. My 
student Scarcelli published a lens-less ghost imaging experiment of chaotic light and raised a 
reasonable question: "Can two-photon correlation of chaotic light be considered as correla- 
tion of intensity fluctuations?" 32J. The lens- less ghost imaging setup of Scarcelli et al. is 
a straightforward modification of the historical Hanbury-Brown and Twiss experiment (HBT) 
33J. Advancing from HBT to the fundamentally interesting and practically useful lens-less 
ghost imaging, what one needs to do is simply move the two HBT photodetectors from far-field 
to near-field. We cannot but stop to ask: What has been preventing this simple move for 50 
years (1956-2006)? Some aspect must be terribly misleading to give us such misled confidence 
not to even try the near-field measurement in half a century. As we know, unlike the first-order 
correlation of radiation that is considered as the interference effect of the electromagnetic waves, 
the second-order correlation of light is treated as statistical correlation of intensity fluctuations. 
Scarcelli et al. pointed out that although the theory of statistical correlation of intensity fluc- 
tuations gives a reasonable explanation to the far-field HBT phenomena, it does not work in 
near-field and consequently does not work for their lens-less near-field ghost imaging experiment 
|34j . It was the idea of statistical correlation of intensity fluctuation that has prevented this 
from happening for 50 years. On the other hand, under the framework of Glauber's theory of 
photodctection, Scarcelli et al. proved a successful interpretation based on the quantum picture 
of two-photon interference. This successes indicates that although the concept of multi-photon 
interference, or the superposition of multi-photon amplitudes, was benefited from the research 
of entangled states, the concept is generally true and applicable to any radiation, including 
"classical" thermal light. Unfortunately, this concept has no counterpart in classical electro- 
magnetic theory of light. Now, we are back to 1996. My student and I have been charged 
with "guilt" again because we have told the physics community a simple truth of the failure of 
a classical idea and adapted the quantum mechanical concept of two-photon superposition to 
"classical" light. 

It was a mistake to keep silence. I have finally resolved to speak about the subject. The 
concept of multi-photon coherence, or the superposition principle of multi- photon amplitudes, 
is important and worthwhile to do, even if I might be burned at the stake. 



Appendix: Presnel propagation-diffraction 

In Fig. \A— 1[ the field is freely propagated from the source plane cto to an arbitrary plane 
a. It is convenient to describe such a propagation in the form of Eq. (|33|) . We now evaluate 
g(K,co; p, z), namely the Green's function for free-space Fresnel propagation-diffraction. 

According to the Huygens-Fresnel principle, the field at a space-time point (p,z,t) is the 
result of a superposition of the spherical secondary wavelets originated from each point on the 
cto plane, see Fig. \A— 1[ 

#>(ftz,t)= f dudKa(uj,K) f dpo^^-e- 1 ^-^" 1 (A-l) 

where A(po) is the complex amplitude, or distribution function, in terms of the transverse 
coordinate po, which may be a constant, a simple aperture function, or a combination of the 
two. In Eq. jA— ID , we have taken zq = and to = on the source plane of <7o as usual. 
In a paraxial approximation, we take the first-order expansion of r' in terms of z and p 

rW* 2 + lp-Pol 2 ^(l + ^y|^). 
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Figure A—l: Schematic of free-space Fresnel propagation. The complex amplitude A(po) is 
composed by a real function A(po) and a phase e~ rK ' Po associated with each of the transverse 
wavevector k on the plane of Co- Notice: only one mode of wavevector k(ft, oS) is shown in 
the figure. 



E( + \p, z, t) is thus approximated as 

#)^«,i)s! / dtodKa(uj,K) f dpo^^-e^e^P-^e-^ 



. ^ . ^ |2 

where e l ~< p ~ Po ' is named as the Fresnel phase factor. 

Assuming the complex amplitude A(po) is composed of a real function A(po) and a phase 
e -iK-p 0j associated with the transverse wavevector and the transverse coordinate on the plane 
of do, which is reasonable for the setup of Fig. — H E(p, z, t) can be written in the following 
form 



E ( -+\p,z,t) = J dujdti a(u,ti)e- lut — J dp A(p )e lR - p0 e l ^ p - M . 
The Green's function u>; p, z) for free-space Fresnel propagation is thus 

g{it,w;p,z) = — [ dpoA(po)e ig -P"G(\p-po\,-). {A-2) 
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In Eq. 2|) we have defined a Gaussian function G(|a|, (3) ~ e 1 ^/ 2 ^" , namely the Fresnel 
phase factor. It is straightforward to find that the Gaussian function G(\a\,f3) has the following 
properties: 

G*(\a\,(3) = G(\a\,-{3), 
G(\d\,p 1 + p 2 ) = G{\d\^x)G{\al(3 2 ), 
G{\a 1 +a 2 \,(3) = G(\d 1 \,P)G(\a 2 \,/3)e i ^ S2 , 

ddG(\a\,/3)e^ s = i^G(\%-~). (A-3) 

P P 

Notice that the last equation in Eq. \A— 3D is the Fourier transform of the G(\a\,f3) function. 
As we shall see in the following, these properties are very useful in simplifying the calculations 
of the Green's functions g(K,uj; p,z). 

Now, we consider inserting an imaginary plane a' between <jq and a. This is equivalent 
having two consecutive Fresnel propagations with a diffraction- free a 1 plane of infinity. Thus, 
the calculation of these consecutive Fresnel propagations should yield the same Green's function 
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as that of the above direct Fresnel propagation shown in Eq. 2| ) : 

g(n,uj;p,z) 

= C 2 — — / dp* dp A(p )G(\p>-pol—)G(\p~p<\,—) 

d\d 2 J a > J ao cdi cd 2 

= c— [ dpoA(p )G(\p-po\,-) (A-4) 

where C is a necessary normalization constant for a valid Eq. J^4 — 4| ) , and z — d\ + di- The 
double integral of dpo and dp' in Eq. jA— 4[ ) can be evaluated as 

/ dp' [ dp A(p )G(\p'-po\,^-)G(\p-p'\,^-) 

Ja' Ja CCii cd 2 

dp A(p )G(p ,^-)G(p,^-) [ dPGip'^ + Ue-^^' 

a c «l c «2 Ja' C "1 «2 

i2vrc d x d 2 f uj > , _ W . ,p p u> d x d 2 

a \ d PoMPo)G(po,-r)G(p,—)G{\— + —\,-( )) 

uj di+ d 2 J ao cdi cd 2 dx d 2 c di + d 2 

i2nc d x d 2 f 2f-\nt\- -I " ^ 



uj di+d 2 J rTQ c(di + d 2 )' 

where we have applied Eq. J A — 3D , and the integral of dp' has been taken to infinity. Substituting 
this result into Eq. 4| ), we thus have 

g{K,u);p,z) = 6 — / dp A(p a ) G{\p - p |, ) 



w di+d 2 Jcro 'c(di+d 2 )' 
= C — [ dp A(p )G(\p-p \,-). 

Z Ja Q CZ 

Therefore, the normalization constant C must take the value of C = — iu> /2uc. The normalized 
Green's function for free-space Fresnel propagation is thus 

g(H, u; p, z) = / dp A(p )G(\p-p \,— ). (-4-5) 

Z7TC 2 /_„ CZ 



References 

[1] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 35, 777 (1935). 
[2] N. Bohr, Phys. Rev. 48, 696, (1935). 

[3] Quantum Theory and Measurement, J. A. Wheeler and W.H. Zurek, eds., Princeton Uni- 
versity Press, Princeton, 1983. 

[4] A. Pais, 'Subtle is the lord... ' The Science and the Life of Albert Einstein, Oxford University 
Press, Oxford and New York, 1982. 

[5] M. DAngelo, Y.H. Kim, S. P. Kulik, and Y.H. Shih, Phys. Rev. Lett. 92, 233601 (2004). 

[6] D. Bohm, Phys. Rev. 85, 166, 180 (1952); D. Bohm, Causality and Chance in Modern 
Physics, D. Van Nostrand Co., Inc., Princeton, 1957; D. Bohm and Y. Aharonov, Phys. 
Rev. 108, 1070 (1957). 



42 



[7] J.S. Bell, Physics, 1, 195 (1964); Speakable and Unspeakable in Quantum Mechanics, Cam- 
bridge University Press, New York, 1987. 

[8] J.F. Clauser and A. Shimony, Rep. Prog. Phys. 41, 1883 (1978). 

[9] A. Aspect, et al., Phys. Rev. Lett. 47 , 460 (1981); A. Aspect, et al, Phys. Rev. Lett. 49, 
91 (1982); A. Aspect, et al, Phys. Rev. Lett. 49, 1804 (1982). 

[10] Y.H. Shih and CO. Alley, Phys. Rev. Lett. 61, 2921, (1988); Z.Y. Ou and L. Mandel, 
Phys. Rev. Lett. 62, 50 (1988); T.E. Kiess, Y.H. Shih, A.V. Sergienko, and CO. Alley, 
Phys. Rev. Lett. 71, 3893 (1993); P.O. Kwiat et al, Phys. Rev. Lett. 75, 4337 (1995). 

[11] E.Schrodinger, Naturwissenschaften 23, 807, 823, 844 (1935); English translations appear 
in ref . [5] . 

[12] D.N. Klyshko, Photon and Nonlinear Optics, Gordon and Breach Science, New York, 1988. 
[13] Y.H. Shih, IEEE J. of Selected Topics in Quantum Electronics, 9, 1455 (2003). 
[14] A. Yariv, Quantum Electronics, John Wiley and Sons, New York, (1989). 
[15] R.J. Glauber, Phys. Rev. 130, 2529 (1963); Phys. Rev. 131, 2766 (1963). 
[16] M. H. Rubin, Phys. Rev. A 54, 5349, (1996). 

[17] J. W. Goodman, Introduction to Fourier Optics, McGraw-Hill Publishing Company, New 
York, NY, 1968. 

[18] Y.H. Shih, IEEE J. of Selected Topics in Quantum Electronics, (2007). 

[19] This effect was first proposed for lithography application, namely quantum lithography, by 
A.N. Boto et al, Phys. Rev. Lett. 85, 2733 (2000). 

[20] M. D'Angelo, M.V. Chekhova, and Y.H. Shih, Phys. Rev. Lett. 87, 013603 (2001). M. 
D'Angelo, M.V. Chekhova, and Y.H. Shih, Phys. Rev. Lett., 87, 013603 (2001). Note: Due 
to the lack of a two-photon absorber, the joint-detection measurement in this experiment 
was on the Fourier transform plane rather than on the image plane. It was implicit in 
Ref. [20] that a second Fourier transform, by inserting a second lens in that experimental 
setup, would transfer the Fourier transform of the object onto its image plane, thus giving 
an image with doubled spatial resolution despite the Rayleigh diffraction limit. It might be 
helpful to point out that the observation of sub- wavelength interference in a Mach Zehnder 
type interferometer cannot lead to sub-diffraction-limited images, except a set of double 
modulated interference pattern. The Fourier transform argument works only for imaging 
setups as is in Ref. [20] 

[21] T.B. Pittman, Y.H. Shih, D.V. Strekalov, and A.V. Sergienko, Phys. Rev. A 52, R3429 
(1995). 

[22] D.N. Klyshko, Usp. Fiz. Nauk, 154, 133 (1988); Sov. Phys. Usp, 31, 74 (1988); Phys. Lett. 
A, 132, 299 (1988). 

[23] M. D'Angelo, A. Valencia, M.H. Rubin, and Y.H. Shih, Phys. Rev. A 72, 013810 (2005). 

[24] D.V. Strekalov, A.V. Sergienko, D.N. Klyshko and Y.H. Shih, Phys. Rev. Lett. 74, 3600 
(1995). 

[25] K.R. Popper, in Open Questions in Quantum Physics, G. Tarozzi and A. van der Merwe, 
eds., D. Reidel Publishing Co., Dordrecht, 1985; K.R. Popper, in Determinism in Physics, 
E.I. Bitsakis and N. Tambakis, eds., Gutenberg Publishing Co., Athens, 1985. 



43 



[26] K.R. Popper, Naturwissenschaften 22, 807 (1934); K.R. Popper, Quan turn Theory and the 
Schism in Physics, Hutchinson, London, 1982. 

[27] For criticisms of Popper's experiment, see for example, D. Bedford and F. Selleri, Lett. 
Nuovo Cimento, 42, 325 (1985); M.J. Collett and R. Loudon, Nature 326, 671 (1987); A. 
Sudbcrg, Philosophy of Science, 52, 470 (1985); A. Sudberg, in Microphysical Reality, A. 
van der Merwe et al, eds., Kluwer Academic, Dordrecht, 1988; M. Home, Experimental 
Metaphysics, R.S. Cohen, M. Home and J. Stachel, eds., Kluwer Academic, Dordrecht, 
1997. 

[28] Y.H. Kim and Y.H. Shih, Foundations of Physics, 29, 1849 (1999). 

[29] D.V. Strekalov, Y.H. Kim, and Y.H. Shih, Phys. Rev. A 60, 2685 (1999). 

[30] C.E. Shannon and W. Weaver, The Mathematical Theory of Communication, University 
of Illinois Press, Urbana, 1949. 

[31] N.J. Cerf and C. Adami, Phys. Rev. Lett. 79, 5194 (1997). 

[32] G. Scarcelli, V. Berardi, and Y.H. Shih, Phys. Rev. Lett., 96, 063602 (2006). 

[33] R. Hanbury Brown and R.Q. Twiss, Nature 177, 28 (1956); 178, 1046, 1447 (1956). 

[34] See "Quantum Imaging" by Y.H. Shih, in sequence with this article. 



44 



